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The  crawling  movement  of  cells  in  response  to  a chemical  gradient  is  a complex 
process  involving  the  orchestration  of  various  intracellular  signaling  molecules.  Although 
a complete  mechanism  for  this  process  remains  elusive,  the  very  first  step  of  gradient 
sensing,  enabling  the  cell  to  perceive  the  direction  of  the  imposed  gradient,  has  become 
more  transparent.  The  increased  understanding  of  this  step  has  been  driven  by  the 
discovery  that  application  of  a weak  chemoattractant  gradient  results  in  the  localization  of 
membrane  phosphoinositides  at  the  front  end  of  the  cell,  which  then  acts  as  a “compass” 
for  the  forward  motility  of  the  cell. 

Cell  migration  plays  a pivotal  role  in  diverse  biological  phenomena  such  as  wound 
healing,  cancer  metastasis  and  inflammatory  response.  Therefore,  an  understanding  of  the 
gradient  sensing  has  vast  applications  in  the  treatment  of  numerous  human  diseases. 

In  this  dissertation,  a mathematical  model  for  gradient  sensing  has  been  formulated. 
This  model  contains  a slow  diffusing  ‘activator’  that  is  synthesized  autocatalytically,  and 
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a fast  diffusing  ‘inhibitor’  that  inhibits  the  activator.  The  activator  and  the  inhibitor  have 
been  identified  from  the  phosphoinositide  signaling  pathways.  Such  a model  explains  the 
localization  of  phosphoinositides  in  response  to  a chemoattractant  gradient  and  also 
simulates  all  the  other  phosphoinositide  dynamics  observed  in  motile  cells. 

The  shape  of  the  steady  state  phosphoinositide  localization  has  been  analyzed  by 
employing  singular  perturbation  theory.  This  technique  exploits  the  difference  in 
diffusion  rates  between  the  activator  and  the  inhibitor  to  analytically  construct  an 
approximation  to  the  phosphoinositide  localization.  Geometric  properties  of  the 
phosphoinositide  localization,  such  as  height  and  width  have  also  been  derived. 

The  stability  of  the  homogenous  phosphoinositide  steady  state  has  been  analyzed. 
This  analysis  shows  that  the  homogenous  steady  state  of  the  phosphoinositides  can 
spontaneously  polarize  into  a phosphoinositide  localization  only  when  the 
phosphoinositide  synthesis  rate  is  at  an  optimum.  This  result  provides  an  experimentally 
testable  characteristic  of  the  model.  Two  types  of  experiments  have  been  done  to  test  this 
hypothesis:  (1)  biophysical  experiments  that  study  the  polarity  at  different  uniform 
chemoattractant  concentrations  and  (2)  biochemical  experiments  that  overexpress  or 
inhibit  an  enzyme  critical  for  phosphoinositide  synthesis.  Results  of  both  these 
experiments  are  consistent  with  model  predictions. 
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CHAPTER  1 

GRADIENT  SENSING  IN  EUKARYOTIC  CELLS 


Introduction 

Directed  cell  migration  accompanies  us  from  conception  to  death. 

• At  birth,  controlled  cell  migration  brings  about  shape  and  organization  to  the 
embryo  [1],  As  the  embryo  develops,  neurons  navigate  to  appropriate  regions  of  the 
body  and  direct  the  formation  of  the  nervous  system. 

• In  an  adult,  cell  migration  is  indispensable  for  the  defense  and  maintenance  of  the 
body.  Neutrophils  rapidly  detect  inflammation  sites  and  engulf  invading  pathogens 
[2].  Fibroblasts  under  the  skin  move  in  large  numbers  towards  wounded  areas  and 
induce  healing  [3]. 

• Cell  migration  can  also  accelerate  death.  Cancer  metastasis  is  caused  by  directed 
cell  migration  of  tumor  cells  from  the  primary  site  to  preferential  sites  of  metastasis 

[4]- 

Hence,  a good  understanding  of  cell  migration  will  be  able  to  shed  light  on  all  these 
biological  phenomena. 


Figure  1-1.  The  five  phases  of  the  chemotactic  cycle  (adapted  from  [5])  (1)  Gradient 

sensing  (2)  Protrusion  of  pseudopod  (3)  Adhesion  to  substratum  (4)  Traction 
of  the  cell  body  and  (5)  Retraction  of  the  cell  tail. 

Most  eukaryotic  cells  move  by  crawling  on  a surface.  This  crawling  movement  is 

initiated  in  response  to  an  external  stimulus,  which  is  frequently  a chemical  concentration 

gradient.  The  resultant  motion  propels  the  cell  forward  along  the  direction  of  highest 
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concentration.  The  chemical  that  induces  the  movement  is  called  chemoattractant  and  the 
movement  itself  is  called  chemotaxis  [6]. 

Eukaryotic  chemotaxis  is  cyclic  and  each  cycle  consists  of  five  distinct  phases 
(Figure  1-1)  [5],  The  cell  first  senses  the  external  chemoattractant  gradient  using 
specialized  receptive  proteins  on  its  membrane.  This  results  in  the  extension  of  a 
protrusion,  called  a pseudopod,  by  polymerizing  actin  in  the  direction  of  maximum 
chemoattractant.  In  order  to  convert  this  response  to  motion,  the  pseudopod  then  adheres 
to  the  substratum.  This  is  then  followed  by  a contraction  phase,  in  which  the  cell  body 
and  nucleus  are  pulled  forward.  The  last  step  in  locomotion  consists  of  two  distinct 
processes,  deadhesion  of  the  protrusion  and  retraction  of  the  tail. 

Each  phase  of  the  chemotactic  cycle  is  a complex  process  involving  the  coordinated 
action  of  a large  constellation  of  signaling  molecules,  many  of  which  have  been 
identified.  Still  lacking,  however,  is  a synthetic  theory  explaining  how  these  molecules 
are  organized  in  space  and  time. 

This  work  is  primarily  concerned  with  understanding  the  first  phase  of  the  cycle, 
gradient  sensing,  the  mechanism  that  enables  the  cell  to  read  the  external  gradient  and 
extend  a pseudopod  precisely  at  its  leading  edge,  the  region  exposed  to  the  highest 
chemoattractant  concentration. 

The  chemoattractant  gradients  imposed  in  the  extracellular  space  are  often  quite 
small  (1-2%  concentration  change  over  the  length  of  the  cell)  [7],  but  the  actin  polymers, 
that  constitute  the  pseudopod,  are  found  exclusively  at  the  leading  edge  [8,  9],  Thus,  a 
key  problem  of  gradient  sensing  is  the  elucidation  of  the  mechanism  that  mediates  the 


3 


formation  of  a highly  polarized  distribution  of  actin  polymers  in  response  to  a relatively 
mild  chemoattractant  gradient. 


Figure  1-2.  Signaling  involved  in  gradient  sensing.  The  chemoattractant  binds  to  G- 
protein  coupled  receptors  on  the  cell  surface,  which  in  turn  induces  the 
dissociation  of  G-proteins.  Active  G-proteins  then  propagates  the  signal 
downstream,  eventually  resulting  in  localized  actin  polymerization  at  the 
leading  edge  of  the  cell.  In  spite  of  a shallow  external  gradient,  the  response 
observed  is  highly  localized  at  the  leading  edge  of  the  cell. 

Several  cell  types  have  been  used  as  model  systems  for  studying  gradient  sensing. 
They  include  fast  moving-cells  such  as  neutrophils  [10]  and  the  slime  mold  amoeba, 
Dictyostelium  [11],  and  slow-moving  cells  such  as  neurons  [12],  budding  yeast  [13,  14] 
and  fibroblasts  [15],  Although  several  aspects  of  this  work  apply  to  slow  moving  cells, 
the  primary  focus  will  be  on  the  fast-moving  cell  types. 

The  chemoattractant  gradient  is  transmitted  to  the  actin  polymerization  machinery 
by  a signal  transduction  pathway  that  starts  with  receptors  on  the  cell  surface  and 
terminates  in  proteins  that  catalyze  actin  polymerization.  In  Dictyostelium  and 
neutrophils  the  receptors  that  regulate  chemotaxis  are  generally  G-protein  coupled 
receptors  (GPCRs).  Upon  ligand  binding,  these  receptors  activate  heterotrimeric  G- 
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proteins  and  dissociate  them  into  Ga  and  GPy  subunits.  Each  of  these  subunits  then 
activates  a train  of  signaling  events,  which  eventually  activates  proteins  that  catalyze 
actin  polymerization. 

It  is,  therefore,  conceivable  that  actin  polymers  inherit  their  highly  polarized 
distribution  from  some  molecule  that  is  upstream  of  the  polymers  in  the  pathway.  Hence, 
it  becomes  crucial  to  identify  the  first  polarized  component  in  the  chemoattractant 
activated  signal  transduction  pathway.  Recent  experiments  [15-18]  have 

• Shown  that  membrane-resident  phosphoinositides,  phosphatidylinositol  3,4,5- 
phosphate  (PI(3,4,5)P3)and  phosphatidylinositol  3,4  phosphate  (PI(3,4)P2)  are 
among  the  earliest  polarized  components  of  the  signal  transduction  pathway 
involved  in  gradient  sensing. 

• Studied  the  spatiotemporal  dynamics  of  these  phospholipids  in  response  to  various 
chemoattractant  profiles. 

Identification  of  the  First  Polarized  Component 

Motile  cells  are  transfected  with  chimeric  proteins  made  by  fusing  a fluorescent 
protein  either  to  the  molecule  of  interest  or  to  a substance  that  binds  specifically  to  the 
molecule  of  interest  and  thus  “reports”  on  it.  For  example,  the  PH  domain  of  Akt  binds 
specifically  to  PI(3,4)P2  and  PI(3,4,5)P3,  so  that  the  fusion  protein  GFP-PH-Akt  reports 
on  the  distribution  of  these  phosphoinositides.  A transfected  cell  is  then  exposed  to  a 
chemoattractant  gradient  by  releasing  chemoattractant  from  a micropipette  and  the 
resultant  intracellular  distribution  of  the  fluorescent  probes  is  visualized  using 
fluorescence  microscopy  (Figure  1-3).  This  approach  allows  the  study  of  spatial  polarity 
at  each  stage  of  the  signal  transduction  pathway  (Figure  1-2). 
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Figure  1-3.  The  principle  of  the  fluorescent  imaging  experiments. 

Specifically,  it  has  been  found  that 

1 Membrane  receptors  are  not  polarized.  It  was  initially  suggested  that  the 

localized  actin  polymerization  could  be  explained  by  a non-uniform  distribution  of 
receptors  on  the  cell  membrane  [19,  20],  However,  recent  studies  have  univocally 
demonstrated  that  when  exposed  to  a chemoattractant  gradient,  receptors  fail  to 
redistribute  and  are  in  fact  uniformly  distributed  [21,  22], 

2.  Receptor  occupancy  and  G-proteins  are  not  significantly  polarized.  Cells 
containing  inactive  G-proteins  fail  to  chemotax  [23],  Hence,  receptors  and  G- 
proteins  are  essential  for  transmission  of  the  chemoattractant  signal.  However, 
when  cells  are  exposed  to  a chemoattractant  gradient,  the  ligand  bound  receptors 
and  active  G-proteins  show  only  a shallow  anterior-posterior  profile  [24,  25],  It 
follows  that  receptors  and  G-proteins  are  required,  but  are  not  the  source  of  the 
localized  amplification  observed  at  the  leading  edge  of  the  cell. 

3.  Phosphoinositides  are  significantly  polarized.  Within  5-10secs  of  applying  a 
chemoattractant  gradient,  PH  domains  of  signaling  molecules  such  as  Akt,  CRAC, 
PhdA  and  Btk  strongly  polarize  at  the  edge  of  the  cell  membrane  that  received  the 
strongest  chemoattractant  stimulus  [16-18],  These  PH  domains  report  specifically 
on  the  intracellular  distribution  of  membrane-resident  phosphoinositides,  PI(3,4)P2 
and  PI(3,4,5)P3.  In  neutrophils,  the  gradient  of  these  PH  domains  is  nearly  six  times 
the  chemoattractant  gradient  [18].  It  follows  that  these  phosphoinositides  are 
among  the  earliest  polarized  components  of  the  signal  transduction  pathway. 
Hereafter,  this  phenomenon  will  be  referred  to  as  phosphoinositide  localization. 

4.  Phosphoinositides  polarize  independent  of  the  actin  cytoskeleton.  The 

phosphoinositide  localization  is  tightly  accompanied  with  actin  polymerization  and 
pseudopod  formation.  This  raised  the  possibility  that  actin  polymerization  is 
necessary  for  the  translocation  of  phosphoinositides.  In  order  to  test  this,  a toxin 
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called  Latranculin  A is  added  to  depolymerize  the  actin  and  give  the  cell  a rounded 
morphology.  However,  cells  still  recruited  the  phosphoinositides  asymmetrically  to 
the  face  closest  to  the  pipette  [11,  16,  18].  This  showed  that  the  localization  of 
phosphoinositides  can  take  place  independent  of  the  actin  cytoskeleton. 

Hence,  the  key  problem  in  the  study  of  gradient  sensing  now  becomes:  What  is  the 

mechanism  underlying  the  phosphoinositide  localization? 


Figure  1-4.  Simple  model  for  localization  of  phosphoinositides.  The  activity  of  the 

enzyme  that  synthesizes  phosphoinositides  is  a sharp  non-linear  function  of 
the  active  receptors.  Any  shallow  chemoattractant  gradient  centered  around 
the  threshold  would  yield  a large  phosphoinositide  gradient.  This  model 
shows  that  the  phosphoinositide  response  is  simply  an  amplified  version  of  the 
chemoattractant  gradient. 

At  first  sight,  the  localization  of  phosphoinositides  seems  explicable  in  terms  of  a 
simple  amplification  model.  It  suffices  to  postulate  that  the  phosphoinositide  synthesis 
responds  to  receptor  activation  in  a highly  cooperative  manner  (Hill-type  kinetics)  (see 
Figure  1-4).  In  this  case,  the  phosphoinositide  distribution  will  be  similar  in  shape,  but 
steeper  in  slope,  when  compared  to  the  chemoattractant  concentration  profile.  In  other 
words,  the  phosphoinositide  distribution  is  an  amplified  version  of  the  chemoattractant 
concentration  profile.  However,  several  experiments  show  that  eukaryotic  gradient 
sensing  is  not  a matter  of  simple  amplification.  In  each  of  the  experiments  below,  the 


distribution  of  the  phosphoinositide  localization  differs  from  the  chemoattractant  profile 
imposed  on  the  cell. 
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Figure  1-5.  Schematic  representation  of  the  spatiotemporal  dynamics.  Upper  panel : 

When  a resting  cell  (shown  in  1)  is  exposed  to  a uniform  increase  in 
chemoattractant  concentration,  it  immediately  responds  by  increasing 
phosphoinositides  all  over  the  membrane  (shown  in  2).  This  is  followed  by  the 
spontaneous  polarization  (shown  in  3)  of  phosphoinositides  along  a random 
direction.  Middle  panel:  When  a resting  cell  (shown  in  4)  is  exposed  to  a 
gradient  of  chemoattractant,  it  accumulates  phosphoinositides  at  the  point  of 
highest  chemoattractant  concentration  (shown  in  5).  When  the  direction  of  the 
gradient  is  switched  the  existing  localization  turns  in  the  direction  of  the  new 
source  of  chemoattractant  (shown  in  6).  This  phenomenon  is  called  polarized 
sensitivity.  Lower  panel:  When  a resting  cell  (shown  in  7)  is  exposed  to  two 
unequal  chemoattractant  sources,  it  responds  to  both  sources  initially. 
However,  as  time  progresses  the  response  to  the  stronger  source  grows  and  the 
response  to  the  weaker  source  is  completely  abolished.  We  refer  to  this  as 
unique  localization. 

Spatiotemporal  Dynamics 

Motile  cells  are  exposed  to  a variety  of  chemoattractant  profiles  which  include  (1) 
steady  or  time  varying  gradients  produced  by  releasing  chemoattractant  from  one  or  more 
chemoattractant  sources  and  (2)  steady  uniform  profiles  obtained  by  immersing  the  cell  in 
chemoattractant.  These  experiments  result  in  the  manifestation  of  the  following 
spatiotemporal  dynamics  (see  Fig  1-5):  (1)  unique  localization  (2)  polarized  sensitivity 
(3)  adaptation  and  (4)  spontaneous  polarization.  In  the  following  section,  these  dynamics 
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will  be  defined  precisely  by  describing  the  experiments  and  the  corresponding 
observations. 

Unique  Localization 

In  normal  motile  cells,  only  one  leading  edge  ultimately  develops,  regardless  of  the 
external  signal.  Foxman  et  al.,  1997,  provided  a convincing  demonstration  of  this 
behavior  by  exposing  a cluster  of  motile  cells  to  two  chemoattractant  sources  located  at 
different  distances  from  the  cluster.  All  the  cells  migrated  towards  the  closer  source, 
showing  that  cells  respond  to  the  stronger  stimulus  and  completely  ignore  the  weaker 
stimulus.  This  phenomenon  is  referred  to  as  unique  localization. 


Figure  1-6.  Dependence  on  pre-existing  polarity.  Response  of  three  cells  labeled  a,  b,  and 
c to  sequential  stimulation  by  two  chemoattractant  sources  (from  [26]).  The 
top  panel  shows  all  three  cells  moving  towards  the  sole  chemoattractant 
source  (pipette  1).  The  bottom  panel  shows  the  response  when  an  additional 
source  (pipette  2)  is  turned  on  some  time  later.  Cell  ‘a’  immediately  turns 
toward  pipette  2,  whereas  cells  ‘b’  and  ‘c’  show  no  response  at  all. 

Polarized  Sensitivity 

If  a cell  that  has  localized  phosphoinositides  in  a certain  direction  is  exposed  to  a 
modest  chemoattractant  gradient  along  a different  direction,  a new  localization  does  not 
develop  at  the  point  with  the  highest  chemoattractant  concentration.  Instead,  the  existing 
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localization  turns  and  reorients  itself  along  the  new  gradient.  This  phenomenon  is  called 
polarized  sensitivity,  since  the  turning  response  suggests  that  pre-existing  leading  edge  or 
“pole”  is  more  sensitive  to  chemotactic  signals  than  all  other  regions  of  the  cell. 
Interestingly,  if  the  new  chemoattractant  gradient  is  relatively  large  and  localized,  the 
existing  pseudopod  retracts  and  a new  one  grows  along  the  direction  of  the  new  gradient. 

The  induction  of  a turning  response  depends  not  only  on  the  external  concentration, 
but  also  on  the  preexisting  polarity  of  the  cell.  This  is  dramatically  illustrated  by  an 
experiment  shown  in  Figure  1-6  [26].  The  upper  panel  shows  three  cells  labeled  a, b and  c 
moving  toward  the  chemoattractant  source  shown  as  pipette  2.  Interestingly,  even  though 
cell  b is  much  closer  to  pipette  2 than  cell  a,  it  does  not  respond  at  all,  whereas  the  more 
distant  cell  a turns  toward  pipette  2.  This  shows  that  the  cells  b and  c that  are  closer  to 
pipette  1,  and  hence  likely  to  be  more  polarized,  are  less  responsive  to  the  subsequent 
influence  of  pipette  2. 

Adaptation  and  Spontaneous  Polarization 

When  resting  cells  are  stimulated  with  a uniform  chemoattractant  concentration,  the 
phosphoinositides  in  the  membrane  increase  uniformly  within  5-10  seconds.  However,  if 
this  uniform  concentration  is  maintained  at  the  same  level,  the  phosphoinositides 
decrease  to  their  prestimulus  steady  state  [11,  16].  This  phenomenon  is  called  adaptation 
[27].  The  mechanism  underlying  adaptation  is  only  partially  understood. 

In  Dictyostelium,  earlier  work  showed  that  adaptation  was  perfect,  i.e.,  the 
membrane  phosphoinositides  returned  to  their  basal  steady  state  [16].  However,  recent 
work  shows  that  cells  frequently  form  multiple  patches  of  phosphoinositide  localization, 
instead  of  completely  adapting  to  their  prestimulus  state,  suggesting  that  adaptation  is  in 
fact  imperfect  [28]. 
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Adapatation  is  manifested  in  a slightly  different  manner  in  neutrophils  [29].  After 
exhibiting  a uniform  initial  increase,  the  phosphoinositides  spontaneously  localize  at  a 
random  location  on  the  membrane  and  decrease  at  all  other  regions  of  the  cell  [18,  30]. 
The  region  of  phosphoinositide  localization  subsequently  becomes  the  leading  edge  of 
the  cell.  This  phenomenon  is  called  spontaneous  polarization  to  emphasize  the  fact  that 
the  phosphoinositides  and  the  cell  morphology  polarize  even  though  the  chemoattractant 
environment  is  macroscopically  uniform  [14].  The  time  taken  for  the  initial,  uniform 
increase  is  independent  of  the  chemoattractant  concentration  [29].  However,  the 
amplitude  of  the  initial  response  and  the  time  taken  for  the  spontaneous  polarization  to 
occur  is  directly  proportional  to  the  chemoattractant  concentration. 

Proposed  model  for  Gradient  Sensing 

From  the  dynamics  discussed  above,  it  is  possible  to  discern  essential  components 
for  a model  of  gradient  sensing.  A large  accumulation  of  membrane  phosphoinositides  at 
the  leading  edge  implies  that  the  external  chemoattractant  signal  is  being  amplified.  The 
lack  of  a response  at  low  concentrations  of  chemoattractant  implies  the  existence  of  a 
threshold  below  which  there  is  no  amplification.  These  properties  are  exhibited  in 
systems  containing  an  activator,  a substance  that  produces  more  of  itself  (ensures 
amplification)  and  possesses  co-operative  kinetics  (ensures  a threshold  for  amplification). 

A model  consisting  of  only  an  activator  cannot  exhibit  a localized  response.  When 
sub-threshold  perturbations  are  applied,  the  system  comes  back  to  the  original 
homogenous  steady  state.  On  the  other  hand  when  supra-threshold,  localized 
perturbations  are  applied  to  such  a system,  the  activator  is  autocatalytically  produced  at 
the  site  of  the  perturbation  and  will  diffuse  freely  in  order  to  ultimately  reach  a uniformly 
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high  activator  profile.  The  localized  response  can  only  be  achieved  by  limiting  the 
uncontrolled  propagation  of  the  activator,  using  a component  called  the  inhibitor. 

This  leads  us  to  consider  a Gierer  and  Meinhardt  type  activator-inhibitor  model 
[31]  that  allows  for  the  formation  of  a stable  pattern.  This  model  consists  of  two  variables 
- an  activator  and  an  inhibitor.  An  activator  enhances  production  of  itself  and  the 
production  of  the  inhibitor.  An  inhibitor  enhances  production  of  itself  and  inhibits  the 
production  of  the  activator.  The  most  general  form  of  the  model  is  as  follows 


da 

d2a 

Activator : 

~dt 

= /(fl,/)  + D.-TT 
ox 

(1.1) 

di 

d2i 

Inhibitor : 

Tt 

-g(a,i)  + Di 

ox 

(1.2) 

where  / (a,  i)  and  g{a,i ) express  the  rate  of  formation  of  the  activator  and  the  inhibitor. 
Da  and  D , denote  the  lateral  diffusivities  of  the  activator  and  inhibitor,  respectively.  Let 

us  assume  that  both  variables  are  initially  at  a homogenous  steady  state  and  that  a supra- 
threshold  perturbation  is  given  at  a random  location  in  space.  If  the  diffusion  coefficients 
of  both  these  components  are  of  the  same  order  ( Da  « D, ),  localized  production  of  the 

activator  is  inhibited  by  the  production  of  inhibitor  at  the  very  same  region.  Again,  a non- 
uniform  steady  state  ceases  to  exist.  Gierer  and  Meinhardt  [3 1]  proposed  that  the 
localization  of  the  activator  is  possible  only  if  the  inhibitor  is  much  faster  diffusing  than 
the  activator  ( D,  » Da ).  In  this  case,  the  inhibitor  would  immediately  diffuse  away  the 

site  of  the  perturbation.  This  rapid  diffusion  of  the  inhibitor  has  a two-fold  effect:  (a)  It 
allows  for  the  amplification  of  the  activator  at  the  site  of  the  perturbation,  (b)  It  increases 
the  inhibitory  effect  away  from  the  perturbation  site,  and  prevents  the  build  up  of  the 
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activator.  This  results  in  the  stable  localization  of  the  activator  at  the  perturbation  site. 
Turing  [32],  was  the  first  to  discover  the  diffusion  driven  instability  of  the  homogenous 
steady  state. 


Table  1-1. 

Models  of  gradient  sensing. 

Spatiotemporal  dynamics 

Model 

Localization 

Adaptation 

Spontaneous 

polarization 

Polarized 

sensitivity 

Unique 

localization 

[33] 

[34] 

S 

V 

[35] 

s 

S 

[36,  37] 

s 

V 

[38] 

s 

V 

[39] 

s 

s 
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The  desire  to  capture  the  spatiotemporal  dynamics  of  the  phosphoinositides  has 
spurred  the  development  of  several  mathematical  models.  Each  of  these  models  present 
elegant  mechanisms  that  can  explain  some  of  the  dynamics  of  gradient  sensing.  However, 
none  of  these  models,  except  for  the  model  described  in  this  work  [39],  have  predicted  all 
of  the  spatiotemporal  dynamics.  Table  1-1  shows  a comparison  of  these  models  with 
respect  to  the  spatiotemporal  dynamics  that  they  exhibit. 

Three  of  these  models  contain  a short-range  activator  that  is  synthesized 
autocatalytically,  and  a long-range  inhibitor  that  inhibits  the  synthesis  of  the  activator 
[33-35],  They  differ  only  with  respect  to  the  postulated  mechanisms  of  the  activation  and 
inhibition  and  the  dynamics  that  have  been  explained.  Two  other  models  contain  a long- 
range  inhibitor  but  no  activator;  i.e.,  there  is  no  autocatalytic  synthesis  [36-38].  These 
two  models  differ  with  respect  to  the  reaction  kinetics  — the  synthesis  rate  of  the 
inhibitor  is  rapid  in  the  first  case,  and  slow  in  the  second  case.  Recently,  Haugh  and 
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coworkers  [40,  41],  have  proposed  a model  for  fibroblast  gradient  sensing  based  on  the 
turnover  and  diffusion  of  a single  intracellular  component. 

This  dissertation  has  been  organized  in  the  following  manner.  In  Chapter  2,  a 
signaling  pathway  called  the  phosphoinositide  cycle  has  been  proposed  for  gradient 
sensing.  Candidates  for  activator  and  inhibitor  have  been  identified  from  this  signaling 
pathway.  Based  on  these  assumptions,  a mathematical  model  for  gradient  sensing  has 
been  formulated.  This  model  has  been  used  to  simulate  the  phosphoinositide  dynamics 
observed  when  the  cell  is  exposed  to  various  types  of  chemoattractant  profiles. 
Specifically,  the  manifestation  of  five  distinct  types  of  dynamics:  (1)  phosphoinositide 
localization  (2)  unique  localization  (3)  polarized  sensitivity  (4)  adaptation  and  (5) 
spontaneous  polarization,  have  been  explained  using  this  model.  These  results  suggest 
that  an  activator-inhibitor  type  mechanism  is  in  fact  consistent  with  the  dynamics  of 
gradient  sensing. 

Chapter  3 describes  the  analytical  construction  of  the  phosphoinositide  localization. 
This  analysis  exploits  the  difference  in  diffusion  rates  of  the  activator  and  the  inhibitor,  to 
construct  a zeroth  order  approximation  to  the  phosphoinositide  localization.  Geometric 
properties  of  the  phosphoinositide  localization,  such  as  height  and  width  have  also  been 
derived.  Chapter  3 also  discusses  the  stability  of  the  homogenous  resting  state  of  the 
phopshoinositides  with  respect  to  small  non-homogenous  perturbations.  This  analysis 
allows  us  to  determine  the  conditions  under  which  the  homogenous  steady  state 
spontaneously  polarizes  into  a phosphoinositide  localization. 

The  results  of  the  model  simulations,  in  Chapter  2,  show  that  the  width  of  the 
phosphoinositide  localization  changes  with  the  rate  of  phosphoinositide  synthesis.  This 
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provides  a simple  way  to  experimentally  test  the  model  that  has  been  formulated  for 
gradient  sensing.  Chapter  4 describes  these  experiments  that  have  been  performed  on  a 
neutrophil  cell  line.  Specifically,  two  types  of  experiments  have  been  performed:  a) 
biophysical  experiments,  where  the  extent  of  polarity  has  been  studied  at  various  mean 
concentrations  of  chemoattractant  and  b)  biochemical  experiments,  in  which  the  activity 
of  the  enzyme  that  synthesizes  phosphoinositides  has  been  manipulated  to  verify  if  the 
rate  of  phosphoinositide  synthesis  can  influence  phosphoinositide  polarization.  The 
results  of  both  these  experiments  turn  out  to  be  consistent  with  the  hypothesized  model. 


CHAPTER  2 

MODEL  FOR  GRADIENT  SENSING 


Many  of  the  signaling  components  involved  in  gradient  sensing  have  been 
identified.  However,  the  exact  mechanism  underlying  this  process  remains  unclear.  The 
simple  model  formulated  in  this  work,  reproduces  nearly  all  the  spatiotemporal  dynamics 
observed  during  gradient  sensing.  The  model  is  based  on  a signal  transduction  pathway 
called  the  phosphoinositide  cycle.  Although  the  role  of  this  pathway  in  gradient  sensing 
has  not  been  comprehensively  established,  the  success  of  this  remarkably  concise  model 
suggests  that  the  molecular  players,  whatever  is  their  identity,  probably  possess  the 
properties  of  the  variables  in  the  model.  Hence,  the  model  can  help  identify  the  type  of 
reaction  network  that  could  give  rise  to  the  observed  dynamics. 

Signal  Transduction  Pathway:  The  Phosphoinositide  Cycle 
The  signaling  pathway  that  follows  receptor  activation  is  the  subject  of  ongoing 
research.  In  Dictyostelium  discoideum  and  neutrophils,  receptor-ligand  binding  activates 
heterotrimeric  G-proteins.  Activated  G-proteins  can  activate  PI3K1  by  direct  binding 
[42],  resulting  in  the  synthesis  of  PIP3  from  PIP2  (see  Figure  2-1).  The  PIP3  thus 
produced  and  the  G[}y  subunit  then  synergistically  activates  the  membrane-resident 
protein,  P-Rexl  [43-45],  which  belongs  to  the  Dbl  family  Rac-GEFs  (guanine-nucleotide 


1 Abbreviations  Used:  PI,  phosphatidylinositol;  PIP,  phosphatidylinositol  4-phosphate;  PIP2, 
phosphatidylinositol  4,5-phosphate;  PIP3,  phosphatidylinositol  3,4,5-phosphate;  DG,  diacylglycerol;  IP3, 
inositol  1,4,5-triphosphate;  PI3K,  phosphatidylinositol  3-kinase;  PI5K,  phosphatidylinositol  4-phosphate  5- 
kinase;  PITP,  phosphatidylinositol  transport  protein;  Prexl,  PIP3-dependent  Rac  exchanger 
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exchange  factors  for  Rac)  that  activate  Rac.  Activated  P-Rexl  then  converts  the  inactive 
Rac-GDP  to  active  Rac-GTP. 


ENDOPLASMIC 

MEMBRANE  CYTOSOL  RETICULUM 


Figure  2-1.  The  phosphoinositide  cycle. 

There  is  growing  evidence  that  Rac  mediates  activation  of  both  PI3K  and  PI5K  [46, 
47],  This  has  led  to  the  suggestion  that  PIP3  and  Rac  function  in  a positive  feedback  loop 
for  more  synthesis  of  PIP2  and  PIP3  [34,  48,  49],  Activation  of  PI5K  creates  yet  another 
positive  feedback  loop  because  this  increases  the  synthesis  rate  of  PIP2  and  its 
downstream  product,  phosphatidic  acid  (PA),  a potent  activator  of  PI5K  [50],  Because  of 
these  two  positive  feedback  loops,  the  synthesis  rate  of  PIP2  and  PfP3  can  rapidly 

accelerate  to  high  levels. 

Such  high  synthesis  rates  of  PIP2  and  PIP3  can  be  sustained  for  no  more  than  a 
second  because  the  concentration  of  phosphatidylinositol  (PI)  in  the  plasma  membrane  is 
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quite  small  [51].  Depletion  of  PI  in  the  plasma  membrane  is  prevented  by  the  cytosolic  PI 
transport  protein  (PITP),  which  transfers  readily  available  PI  from  the  endoplasmic 
reticulum  to  the  plasma  membrane  [52].The  PIP2  formed  by  successive  phosphorylation 
of  PI  is  hydrolyzed  by  PLC  to  diacylglycerol  (DG)  and  cytosolic  inositol  1,4,5- 
triphosphate  ( IP3).  Diacylglycerol  is  converted  to  PA  and  transferred  to  the  endoplasmic 

reticulum  for  regeneration  of  PI.  Inositol  produced  by  rapid  dephosphorylation  of  IP3  via 
multiple  pathways  [53],  also  participates  in  PI  regeneration. 


Figure  2-2.  The  kinetic  scheme  of  the  model.  The  receptors  receive  the  chemoattractant 
signal  and  instantly  transmit  it  to  the  adaptation  subsytem  consisting  of  two 
hypothetical  variables  U and  V . The  variable  U instantly  transmits  the 
signal  to  the  polarization  subsystem,  which  is  an  abstraction  of  the 
phosphoinositide  cycle.  P and  Ps  denote  the  phosphoinositide  pools  in  the 

plasma  membrane  and  the  endoplasmic  reticulum,  and  / denotes  the  inositol 
phosphates  in  the  cytosol.  These  rapid  events  are  followed  by  the  slow  growth 
of  V which  inhibits  the  synthesis  of  U and  restores  it  to  to  its  pre-stimulus 
value. 

Several  observations  suggest  a causal  link  between  PIP2  [54-56]/PIP3  [16-18] 
synthesis  and  lamellipod  extension.  Lamellipods  are  formed  precisely  at  the  same 
location  as  PIP2  [54]/ PIP,  [16,  18]  production.  The  enzymes  that  catalyze  the  synthesis  of 
these  phosphoinostides,  PI3K  [57]  and  PI5K  [58],  are  also  recruited  to  the  leading  edge 
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of  the  cell.  In  addition,  it  has  been  shown  that  PIP2 , in  conjunction  with  GTP-bound 

Cdc42,  is  a strong  activator  of  N-WASP,  which  in  turn  activates  Arp2/3  [59],  Similarly, 
Rac  activation  and  PIP3  production  activate  WAVE  proteins,  which  can  also  activate 
Arp2/3  [59].  Activated  Arp2/3  mediates  actin  polymerization  by  nucleating  the  sides  of 
pre-existing  actin  filaments  [60].  Actin  polymerization  by  Arp2/3  is  believed  to  drive 
lamellipod  protrusion  [61].  Taken  together,  these  facts  suggest  that  the  localization  of 
PIP2  and  PIP,  resulting  from  the  gradient  sensing  mechanism  plays  a crucial  role  in  the 
subsequent  extension  of  the  lamellipod. 

Model 

Figure  2-2  shows  the  structure  of  the  model.  The  extracellular  signal  is  received  by 
the  receptor  and  transmitted  sequentially  to  the  adaptation  and  polarization  subsystems. 

In  earlier  work  [34],  adaptation  was  modeled  by  assuming  slow  desensitization  and 
resensitization  of  receptors  In  view  of  the  subsequent  observations  concerning  the  locus 
of  adaptation  in  the  signaling  pathway  [25,  62],  it  has  been  assumed  that  the  adaptation 
subsystem  lies  below  the  receptors  and  above  the  polarization  subsystem. 
Receptor-ligand  Binding 

It  is  assumed  that  diffusion  of  receptors  is  negligibly  small  and  receptor-ligand 
binding  follows  the  mass  action  law.  The  concentration  of  ligand-bound  (active) 
receptors,  denoted,  is  then  given  by 

— = k+(r-r)l-k~r  (2.1) 

dt  ' 

where  k+  and  k~  are  the  rate  constants  for  association  and  dissociation,  respectively;  rt 
is  the  total  concentration  of  receptors;  and  l is  the  ligand  concentration. 
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Adaptation  Subsytem 

The  molecular  mechanism  of  adaptation  is  unknown.  However,  since  the  interest 
lies  in  studying  the  effect  of  adaptation  on  the  dynamics  of  amplification,  a simple  ‘toy’ 
adaptation  model  formulated  by  [27]  will  suffice.  This  model  encapsulates  the  key 
properties  of  adaptation.  It  consists  of  two  variables  p and  3 that  do  not  diffuse  and 
evolve  according  to  the  equations 


where  p,  3 are  the  dimensionless  concentrations  of  two  hypothetical  entities,  U,  V ; 
p = rtr,  is  the  fraction  of  ligand-bound  receptors;  1/a,  1/6  are  the  time  constants  for  the 
evolution  of  p and  3 . It  is  assumed  that  the  dynamics  of  U are  much  faster  than  the 
dynamics  of  V , i.e.,  a»ft.  Equations  (2.2-2. 3)  comprise  the  adaptation  subsystem  of 
the  model.  This  subsytem  is  only  capable  of  predicting  perfect  adapation.  Although  the 
adaptation  subsystem  has  no  biochemical  relevance,  it  can  effectively  simulate  the 
dynamics  of  adaptation.  Given  any  active  receptor  distribution,  p , the  steady  state 
concentration  of  U always  approaches  zero.  Inasmuch  as  the  polarization  subsystem 
“sees”  the  extracellular  environment  through  the  variable  U (Figure  2-2),  it  becomes 
blind  to  it  after  a sufficiently  large  time  interval.  The  effects  of  such  a property  on  the 
dynamics  of  polarization  will  be  studied. 

Polarization  Subsytem 

To  capture  the  spatiotemporal  dynamics  described  earlier,  the  adaptation  subsystem 
will  be  coupled  to  the  polarization  subsystem.  The  polarization  subsystem  is  an 


(2.2) 


(2.3) 
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abstraction  of  the  phosphoinositide  cycle  (Figure  2-2).  It  contains  three  variables 
corresponding  to  three  “lumped”  pools,  namely,  membrane  phosphoinositides  ( P ), 
cytosolic  inositol  and  its  phosphates  (I),  and  phosphoinositides  in  the  endoplasmic 
reticulum  ( Ps ).  The  concentrations  of  these  variables  are  denoted  by  p , i and  ps , 
respectively.  It  is  assumed  that 

• The  cell  is  two-dimensional  and  disk-shaped.  Thus,  p and  ps  are  based  on  the 
length  of  the  plasma  membrane  and  / is  based  on  the  area  of  the  cytosol. 

• Radial  gradients  of  the  cytosolic  inositol  phosphate  pool  are  negligibly  small.  This 
is  reasonable  because  inositol  phosphates  diffuse  rapidly.  However,  angular 
gradients  of  this  pool  are  considered  because  membrane  phosphoinositides,  being 
slow  diffusers,  develop  steep  angular  gradients.  These  sharp  gradients  could  induce 
mild  angular  gradients  of  inositol  phosphates. 

• There  is  basal  synthesis  and  degradation  of  P and  / . Basal  synthesis  rates  of  P 
and  I , denoted  cp  and  c, , follow  zero-order  kinetics.  Basal  degradation  rates, 

denoted  rpd  and  rid,  obey  first-order  kinetics  with  rate  constants  kp  and  kt, 
respectively,  i.e.. 


rP,d=kPP,^d^k,i 

The  receptor-mediated  rate  of  formation  of  membrane  phosphoinositides  per  unit 
length  of  membrane  is  given  by 


rpj  =kfPP2Ps 

where  kf  is  the  rate  constant  and  p is  the  dimensionless  concentration  of  U . Here, 
we  have  assumed  that  U induces  phosphoinositide  synthesis  by  activating  PI3K. 
The  dependence  on  p2  represents  the  autocatalytic  and  cooperative  kinetics  with 
respect  to  membrane  phosphoinositides.  This  is  an  idealization  of  the  two  positive 
feedback  loops  shown  in  Figure  2-1.  The  membrane  phosphoinositides  will  play  the 
role  of  a local  activator  in  the  model. 

The  inositol  phosphate  pool  ( / ) stimulates  transfer  of  phosphoinositides  from  the 
plasma  membrane  to  the  endoplasmic  reticulum.  The  rate  of  removal  of  membrane 
phosphoinositides  per  unit  length  of  the  membrane  is  assumed  to  be 


rP,r  = Kp * 
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where  kr  denotes  the  rate  constant.  The  rationale  for  this  assumption  is  as  follows. 
Inositol  reacts  with  CDP.DG  to  regenerate  PI  in  the  endoplasmic  reticulum  (Figure 
2-1).  An  increase  in  inositol  levels,  therefore,  promotes  this  reaction  and  drives  the 
transfer  of  PA  from  the  plasma  membrane  to  the  endoplasmic  reticulum.  Thus,  the 
rate  of  removal  of  membrane  phosphoinositides  depends  on  p and  should  also 
depend  on  the  concentration  of  the  inositol  phosphates,  i . The  inositol  phosphates 
will  act  as  a global  inhibitor. 

The  dynamics  of  P , Ps  and  I are  then  governed  by  the  equations 


dp  2 . D d2p 

— -kfup  p-krpi  + cn-knp-\ — r — — 
dt  Vs  rF  p Py  r2  d6l 

(2.4) 

dPs  = (kfpp2p  k pi + c k p)+  Dp‘  8 P° 
dt  ' f p p ' R2  dO2 

(2.5) 

— = s(kfpp2p  -kpi)  + ct- k,i  + ‘ 

dt  V f s r ’ ' ' R2  dO2 

(2.6) 

where  Dp , Dp  and  Z),  denote  the  lateral  diffusivities  of  P , Ps  and  I , respectively,  and 


R denotes  the  cell  radius.  The  factor  s , which  denotes  the  membrane  length  per  unit  cell 
area,  is  required  since  synthesis  and  removal  rates  of  P are  based  on  the  length  of  the 
plasma  membrane. 

The  complete  model  accounting  for  receptor-ligand  binding,  adaptation  and  spatial 
sensing  consists  of  equations  (2. 1-2.6).  It  remains  to  specify  the  boundary  and  initial 
conditions. 

Boundary  and  Initial  Conditions 

Since  the  cell  is  disc-shaped,  all  the  variables  satisfy  periodic  boundary  conditions;  that  is 
where  x = r,p,&,p,ps,i . For  future  use,  note  that  (2.4-2. 5)  and  the  boundary  conditions 


imply  that  the  total  amount  of  phosphoinositides  in  the  plasma  membrane  and  the 


(•2x3.1416 

endoplasmic  reticulum,  I (p  + ps)RdO , is  conserved.  Hence,  the  average 
concentration  of  phosphoinositides  in  the  cell 

1 (-2x3.1416 

p,  = — I (p  + ps)Rd6,  C = circumference  of  the  cell 

C * 

is  constant.  The  initial  distribution  of  a variable  will  be  denoted  by  appending  the 
superscript  ‘-’to  the  variable.  Thus 


r(O,0)  = r ,p(Q,0)  = /U,i9(0,#)  = 3~ 
p(0fl)=p-,ps  (O,O)=pjj(O,0)=i-  ,0<#<2x3. 1416 


Various  initial  distributions  have  been  used  in  the  simulations.  They  are  specified  before 
the  description  of  each  of  the  simulations. 

Dimensionless  Equations 

Using  the  dimensionless  variables 


a „sp_  n =Ps_  lsj_  6 
P r,  ’ P,’  5 P,*  sp,  2x3.1416 


,r  = 


\!{krspt) 


and  dimensionless  parameters, 
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the  following  dimensionless  equations 
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di  _ 2 d2i 


dimensionless  boundary  conditions 


x(0,r)  = x(l,r); 


ck(0,r)  _ c!x(l,r) 


^ 3<f 


x = p,p,u,n,ns,i 


and  dimensionless  initial  conditions  can  be  obtained. 


P(0,  £)  = P~,  M(  0,  £)  = ff,  *9(0,  £)  = 9~ , n{  0, £)  = n~ 


Model  Equations 

The  activation  of  receptors  and  PI3K  is  assumed  to  fast  compared  to  amplification, 
i.e.,  Tr,tu  «:  1 . Hence,  p and  p instantly  achieve  the  quasisteady  state  concentrations. 

The  dynamics  for  all  but  the  smallest  initial  time  intervals  are  then  approximated  by  the 
equations 


X 


(2.7) 


(2.8) 


(2.9) 


— ~Kfun  n-m  + u/„-K„n  + d„ — - 
dr  f s Yp  p p d£2 


(2.10) 
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^s-  = ~(Kfpn2n  -m  + y -k  7r)  + S 

Qt  V S T p p ) p,  q^2 


di  2 d i 

— = k fun  n-m  + u/,- k,i  + o,  — r 
dr  f s ' 3#2 


the  boundary  conditions 


,A  x ,,  . Sx(0,r)  3x(l  ,r) 

*(0,r)  = x(l,r), — — = — ,t  > 0,x  = o,n,ns,i 

o # oq 


and  initial  conditions 


(2.11) 

(2.12) 


(2.13) 


.9(0,#)  = = n~,ns(U)  = <,K<U)  = <~,0  < # < 1 (2.14) 


Simulations 

The  simulations  have  been  done  using  the  NAG  subroutine  D03PHF  [63].  The 
parameter  values  used  in  the  simulations  are  shown  in  Table  2-1.  The  rationale  for  the 
choice  of  parameter  values  can  be  found  in  Appendix  B.  To  facilitate  comparison  of  the 


simulations  with  experimentally  observed  dynamics,  it  is  useful  to  note  that  krsp , ~ l/sec. 


Table  2-1.  Dimensionless  parameter  values  used  in  simulations. 


Tr  = 0. 1 

© 

II 

St 

rv  = 1000 

*,=1 

^ =0.001 

sPt  =0.001 

Kf  = 6 

K,  =0.1 

(/A  =0.01 

fa  = 7.67x1  O'7 

*,=0.1 

o 

© 

II 

Therefore,  each  unit3  of  the  dimensionless  time,  r , corresponds  to  roughly  lsec. 
Since  the  time  scale  of  adaptation  is  on  the  order  of  minutes,  we  assume  that  r , i.e., 
adaptation  is  slow  compared  to  the  reactions  involved  in  amplification. 


2 The  choice  of  parameter  values  and  the  nomenclature  for  the  greek  symbols  are  discussed  in  the 

Appendix. 
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Hence,  the  motion  after  PI3K  activation,  which  is  governed  by  equations  (2.7- 
2.14),  can  be  decomposed  into  two  phases. 

1 . During  the  first  phase,  there  is  no  adaptation.  Thus,  9 remains  constant  at  its 
initial  value,  9~ , and  n,ns,i  evolve  towards  their  quasisteady  state. 

2.  During  the  second  phase,  the  effects  of  adaptation  become  perceptible,  so  that  9 
moves  slowly  towards  its  steady  state  value  p , while  n,ns,i  remain  at  their 
quasisteady  state  values. 

These  shall  be  referred  to  as  the  fast  response  and  slow  response,  respectively.  Five 
distinct  types  of  spatiotemporal  dynamics  were  described  earlier.  Three  of  them, 
phosphoinositide  localization,  polarized  sensitivity  and  unique  localization,  can  be 
completely  understood  by  confining  attention  to  the  fast  response.  To  understand  the 
remaining  two  responses,  adaptation  and  spontaneous  polarization,  the  entire  transient, 
consisting  of  the  fast  and  the  slow  response,  must  be  considered. 

Fast  Responses 

The  fast  response  is  described  by  the  equations 


A Q- 

M = , . $ 

1 4*  X 

(2.15) 

dn  2 _ d27t 

— -Kfpn  n -m  + u/„ -K„n  + d„ — - 

dr  f p p p d£2 

(2.16) 

dii'  i 2 \ - d1n. 

d*=  ( Kfpnns  m + y,p  Kpn)  + 8Pi  ^ 

(2.17) 

di  2 * 

— = Kfun  n-ni  + w.-K.i  + o, — 7 
dr  f s nil 

(2.18) 

with  boundary  conditions 


3 In  Dictyostelium,  adaptation  is  observed  within  50secs[16].  This  time  scale  can  be  accommodated  by 
assuming  krsp,  ~ 10 /sec.  In  this  case,  each  unit  of  r corresponds  to  roughly  O.lsec. 
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x(0,r)  = *(l,r), 


9x(0,r) 


dx(\,T ) 

,T>0,x  = n,7ts,i 


and  initial  conditions 


*(0.£)  = ^",^(0,#)  = n~,i( 0,£)  = r,o  < 4 < 1 


(2.19) 


(2.20) 


Figure  2-3.  Phosphoinositide  localization. (a)  Chemoattractant  concentration  profile  at 
t > 0 (b)  Development  of  a pronounced  membrane  phosphoinositide  peak  at 
4 = 1/2  which  corresponds  to  the  leading  edge  of  the  cell  (c)  Growth  of  the 
uniformly  distributed  inositol  phosphate  pool  (d)  Initial  dynamics  of 
membrane  phosphoinositides  at  any  point  of  the  cell  membrane  with  fi  at  that 
point  as  a control  parameter.  When  n is  small,  the  existence  of  a stable  lower 
steady  state  prevents  amplification  of  all  sufficiently  small  perturbations  of 
/r . At  sufficiently  large  values  of  n , the  lower  and  intermediate  steady  states 
vanish  resulting  in  the  disappearance  of  the  threshold  and  the  amplification  of 
M- 
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Phosphoinositide  localization 

The  experiments  show  that  when  a resting  cell  is  exposed  to  a chemoattractant  gradient, 
phosphoinositides  polarize  within  a few  seconds.  To  simulate  this  experiment,  assume 
that  at  r < 0 , the  cell  is  at  the  resting  steady  state,  ie.,  the  steady  state  obtained  in  the 
absence  of  exogenous  chemoattractant  (A  = 0 ).  Hence,  ju~,  3~  = 0 and  are 

constants  satisfying  the  equations 

0 = i u„-ni-  Kn  n 

0 — - m-K'i 

7TS  =1  -It 

At  r > 0,  a chemoattractant  gradient,  2,(£) , is  applied.  The  subsequent  evolution  of  n , 
ns  and  i is  then  obtained  by  integrating  (2.15-2.20)  with  3~  ,n~  ,n~  ,f  as  determined  by 
the  foregoing  equations. 

Figure  2-3a  shows  the  chemoattractant  profile,  2.(£) , imposed  at  r > 0 . Despite  the 
mild  chemoattractant  gradient,  the  simulation  shows  that  a pronounced  phosphoinositide 
peak  develops  at  £ = 1/2  which  corresponds  to  the  leading  edge  of  cell  (Figure  2-3b). 
Compared  to  the  strongly  polarized  distribution  of  membrane  phosphoinositides,  the 
concentration  profile  of  inositol  phosphates  is  virtually  flat  (Figure  2-3c). 

An  intuitive  explanation  of  these  dynamics  can  be  given  by  observing  that  when  the 
cell  is  at  its  resting  steady  state,  there  is  a threshold  for  peak  formation.  To  see  this, 
perturb  the  profile  of  n from  its  initial  uniform  distribution,  , to  a nonuniform 
distribution,  . Immediately  after  the  perturbation,  ns  - n~  = 1 - n and  i = f . 
Moreover,  since  the  initial  concentration  gradients  are  small,  diffusion  plays  no  role. 
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Hence,  the  initial  phosphoinositide  dynamics  at  any  point  of  the  membrane  may  be 
approximated  by 

dn  , 

— = KfHn  {\-n)-nr + i/yp-Kpn  (2.21) 

Figure  2-3d  shows  the  dynamics  of  this  equation  for  various  values  of  /u . If  /r  is 
small,  there  is  either  a low  steady  state  or  bistable  steady  states.  In  the  first  case,  n is 
stably  maintained  at  the  lower  steady  state  and  in  the  second  case,  the  intermediate  steady 
state  acts  as  a threshold  that  prevents  a significant  increase  in  n . However,  if  fj.  is  large 

enough  at  a point,  the  two  lower  steady  states  disappear  and  every  n~  gets  amplified  to 
the  upper  steady  state.  Similarly,  an  increase  in  i decreases  n to  the  lower  steady  state 
and  a decrease  in  i allows  for  the  amplification  of  n to  the  upper  steady  state. 

The  formation  of  the  phosphoinositide  peak  can  now  be  explained  as  follows.  If 
/*(£)  is  sufficiently  large,  the  threshold  vanishes  and  the  concentration  of  P starts 
increasing  at  all  points  on  the  membrane.  However,  this  increase  in  P does  not  persist  at 
the  trailing  edge  of  the  cell.  It  is  contained  by  the  rapid  diffusion  of  the  cytosolic  inositol 
phosphates  (I)  being  produced  at  the  leading  edge.  The  rapid  diffusion  of  I has  a two- 
fold effect.  In  the  neighborhood  of  the  leading  edge,  the  inhibitory  effect  of  I is 
diminished.  In  other  words,  I is  not  high  enough  to  inhibit  the  autocatalytic  formation  of 
P at  the  leading  edge.  Outside  this  neighborhood,  the  higher  concentration  of  / 
promotes  transfer  of  membrane  phosphoinositides  from  the  plasma  membrane  to  the 
endoplasmic  reticulum.  The  net  effect  of  this  transfer  is  to  deplete  the  plasma  membrane 
of  its  phosphoinositides,  thus  preventing  the  peak  from  spreading  beyond  the  leading 
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edge.  Hence,  within  the  leading  edge,  the  steady  state  concentration  of  P is  higher  than 


the  initial  basal  level.  Outside  the  leading  edge,  it  is  lower  than  the  initial  basal  level. 


Figure  2-4.  Unique  localization  (a)  Chemoattractant  profile  at  r > 0 . The  source  S2  is 
slightly  stronger  than  the  source  . (b)  The  evolution  of  the  corresponding 
phosphoinositide  distribution.  A pronounced  phosphoinositide  peak  develops 
at  the  site  of  source  S2 . 

Unique  localization 

The  experiments  show  that  cells  can  discriminate  between  two  chemoattractant 
sources.  Specifically,  when  cells  encounter  two  chemoattractant  sources  of  unequal 
strengths,  they  migrate  toward  the  stronger  source.  The  simulations  show  that  the  model 
cell  discriminates  between  two  sources  of  almost  identical  strengths. 

To  this  end,  assume  that  a resting  cell  is  exposed  to  a chemoattractant  profile, 

/l(£) , arising  from  two  chemoattractant  sources,  S{  and  S2 , of  almost  identical  strengths 


located  at  and  , respectively  (Figure  2-4a).  The  simulation  shows  an  initial  increase 


in  membrane  phosphoinositides  at  both  sites  of  stimulation  (Figure  2-4b).  However,  as 
time  progresses  only  one  phosphoinositide  peak  persists,  and  it  is  at  the  site  of  the 
stronger  source,  S2 . At  the  site  of  the  weaker  source,  S{ , the  initial  buildup  of 
phosphoinositides  is  completely  extinguished. 
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The  practical  implication  of  this  result  is  the  following.  In  reality,  it  is  unlikely  that 
a cell  will  encounter  two  perfectly  identical  chemoattractant  sources.  The  model  therefore 
implies  that  regardless  of  the  chemoattractant  stimulus,  the  cell  decides  to  migrate  in  one 
and  only  one  direction. 

In  terms  of  the  model,  this  property  of  unique  localization  can  be  explained  as 
follows.  Exposure  to  two  unequal  sources  results  in  an  initial  surge  of  phosphoinositide 
synthesis  at  both  sites  of  stimulation.  However,  since  the  source  at  S2  is  stronger, 

phosphoinositide  synthesis  at  S2  is  larger  than  phosphoinositide  synthesis  at  Sx . At  the 
same  time,  there  is  more  synthesis  and  diffusion  of  the  inositol  phosphates  ( I ) from 
towards  . The  net  effect  of  I ’s  fast  diffusion  towards  is  to  increase  the  threshold  at 
and  thus  prevent  peak  formation. 

This  model  exhibits  unique  localization  for  the  parameter  values  chosen  in  Table  2- 
1.  However,  this  property  depends  on  the  parameter  values  associated  with  the  reaction 
and  diffusion  of  the  activator  and  the  inhibitor.  Variations  in  these  parameter  values  can 
yield  multiple  localizations. 

The  quasisteady  state  as  a function  of  the  chemoattractant  stimulus 

More  insight  into  the  model  can  be  gained  by  determining  the  properties  of  the 
quasisteady  states  corresponding  to  various  chemoattractant  profiles.  Figure  2-5a  shows 
three  chemoattractant  profiles  with  the  same  mean  value  ( Am  = 0.26 ) but  very  different 
slopes.  Figure  2-5b  shows  the  bifurcation  curves  corresponding  to  these  chemoattractant 
profiles.  The  bifurcation  curves  were  obtained  by  varying  the  mean  value,  Am , of  the 
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chemoattractant  profiles.  Thus,  they  reflect  the  quasisteady  states  obtained  when  the 


chemoattractant  profiles  are  “lifted”  vertically  without  undergoing  any  change  of  shape. 


Figure  2-5.  Variation  of  the  quasi-steady  state  as  a function  of  chemoattractant  stimulus. 
(a)Three  chemoattractant  profiles  with  the  same  mean  value  but  different 
gradients.  (b)Variation  of  the  maximum  value  of  the  corresponding 
quasisteady  state  phosphoinositide  distributions  with  respect  to  the  mean 
value  of  the  chemoattractant  profiles.  (c)The  quasisteady  state 
phosphoinositide  distributions  corresponding  to  three  different  mean  values.  If 
the  mean  value  lies  in  Regions  A or  C,  the  phosphoinositide  distribution  is 
nearly  flat.  If  the  mean  value  lies  in  Region  B,  the  phosphoinositide 
distribution  is  strongly  polarized. 


We  begin  by  examining  the  bifurcation  curve  corresponding  to  the  uniform 
chemoattractant  distribution  (Figure  2-5 b,  dotted  curve).  Such  uniform  distributions 
result  in  homogeneous  quasisteady  states  that  satisfy  equations  (2.15-2.20)  with  k = km , 
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= 0 and  8p  = 8p  = 8(=0.  The  quasisteady  state  phosphoinositide  concentration 
increases  gradually,  and  then  saturates  at  sufficiently  large  values  of  Xm . 

Next,  consider  the  quasisteady  states  corresponding  to  the  nonuniform 
chemoattractant  distributions.  The  curves  in  Figure  2-5b  and  c show  that  when  the  mean 
values  of  the  nonuniform  chemoattractant  distributions  are  small  or  large,  both  the 
maximum  and  the  distribution  of  the  quasisteady  state  phosphoinositide  concentration 
become  very  similar  to  those  observed  in  response  to  uniform  chemoattractant 
distributions.  However,  there  is  an  intermediate  range  of  mean  values  over  which  the 
maxima  and  the  distribution  of  the  phosphoinositide  concentration  diverge  dramatically 
from  their  homogeneous  counterparts.  For  mean  values  lying  within  this  intermediate 
range,  the  nonhomogenous  quasisteady  state  phosphoinositide  distributions  are  highly 
nonuniform. 

These  results  can  be  rationalized  as  follows.  As  A increases,  so  does  the 
parameter,  tcf(i , which  is  the  ratio  of  the  characteristic  velocities  of  phosphoinositide 
synthesis  and  removal.  If  Xm,  and  hence,  Kffj  is  small,  the  inhibitory  action  of  the 
inositol  phosphates  is  so  fast  compared  to  phosphoinositide  synthesis  that  the 
phosphoinositide  peak  is  abolished  before  it  can  form.  If  Xm  is  large,  the  inositol 

phosphate  pool  responds  so  slowly  compared  to  phosphoinositide  synthesis  that  the  peak 
spreads  across  the  entire  cell  before  the  inhibition  can  contain  it.  In  either  case,  the 
quasisteady  state  phosphoinositide  distribution  is  nearly  uniform.  It  is  only  at 
intermediate  values  of  Am  that  the  time  scales  of  phosphoinositide  synthesis  and 
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inhibition  are  comparable,  and  nonuniform  phosphoinositide  distributions  can  be 
sustained. 

Polarized  sensitivity 

The  model  has  also  been  used  to  explore  the  phenomenon  of  polarized  sensitivity. 
To  this  end,  a resting  cell  is  allowed  to  form  a quasisteady  state  phosphoinositide  peak  in 
response  to  a chemoattractant  gradient,  A,(£) . This  quasisteady  state  is  the  initial 

condition  of  the  cell.  Thus,  >9~=0  and  satisfy  the  boundary  value  problem 

defined  by  the  equations 


and  the  boundary  conditions  (2.19).  After  the  quasi-steady  state  is  reached,  the  polarized 
cell  is  subjected  to  a new  chemoattractant  gradient,  * A,(£)  ■ The  subsequent 

evolution  of  P , Ps  and  I in  response  to  this  new  gradient  is  obtained  by  integrating 
(2.15-2.20)  with  A = X1(%). 

The  simulations  show  that  if  the  new  chemoattractant  gradient  declines  gradually 
from  its  maximum  in  such  a way  that  the  receptors  in  the  neighborhood  of  the  pre- 
existing phosphoinositide  peak  sense  the  influence  of  the  new  gradient,  the  pre-existing 
peak  moves  like  a traveling  wave  to  the  point  at  which  the  new  chemoattractant  gradient 
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Figure  2-6.  Response  of  a pre-existing  phosphoinositide  peak  to  a new  gradient,  (a)  The 
active  receptor  concentration  (b)  When  the  new  gradient  is  shallow,  the 
phosphoinositide  peak  moves  as  a traveling  wave  towards  the  new  steady  state 
(c)  If  the  new  gradient  is  localized  and  large  enough,  the  original 
phosphoinositide  peak  collapses  and  a new  peak  is  created  at  the  maximum  of 
the  new  gradient. 

has  a maximum  (Figure  2-6b).  On  the  other  hand,  if  the  new  chemoattractant  profile  is 
large  enough  to  overcome  the  threshold  and  is  localized  such  that  the  pre-existing  leading 
edge  does  not  sense  the  new  gradient,  then  there  is  no  wave  motion.  Instead,  the  pre- 
existing peak  retracts,  and  a new  peak  grows  at  the  maximum  of  the  new  gradient  (Figure 
2-6c).  These  results  are  consistent  with  the  phosphoinositide  dynamics  observed  in 
response  to  changes  in  chemoattractant  gradients. 
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To  explain  the  wave-like  motion  of  the  peak  shown  in  Figure  2-6b,  it  is  useful  to 
observe  that  the  quasi-steady  phosphoinositide  peak  formed  in  response  to  the  first 
gradient  is  “inert”  everywhere  except  in  the  two  thin  “transition  layers”  surrounding  the 
peak  within  which  there  is  a sharp  change  in  the  gradient  of  membrane 
phosphoinositides.  By“Inert”,  it  means  that  outside  these  transition  layers,  nothing  is 
happening  at  steady  state  - there  is  neither  diffusion  nor  synthesis  of  membrane 
phosphoinositides.  The  transition  layers,  on  the  other  hand,  are  sites  of  intense  activity.  In 
the  upper  half  of  a transition  layers,  there  is  rapid  synthesis  of  membrane 
phosphoinositides  which  then  diffuse  into  the  lower  half  of  the  transition  layer,  from 
where  they  are  promptly  removed.  The  quasi-steady  state  is  maintained  by  this  precarious 
balance  between  synthesis  of  membrane  phosphoinositides  in  the  upper  half  of  the 
transition  layers  and  their  removal  in  the  lower  half  of  the  transition  layers.  If  the  balance 
is  disturbed  by  imposing  a shallow  chemoattractant  gradient  that  increases  the  rate  of 
phosphoinositide  synthesis  relative  to  its  rate  of  diffusion,  the  transition  layer  moves  in  a 
wave-like  fashion  at  a velocity  that  is  proportional  to  the  net  rate  of  phosphoinositide 
accumulation  within  the  transition  layer. 

The  response  to  steep  chemoattractant  gradients  shown  in  Figure  2-6c  can  be 
explained  as  follows.  After  the  quasi-steady  state  has  developed,  the  concentration  of  the 
inhibitor,  I , is  high  throughout  the  cell.  This  tends  to  increase  the  threshold  at  all  points 
of  the  plasma  membrane.  When  the  gradient  is  switched,  the  active  receptor 
concentration  decreases  at  the  previous  “front”  and  increases  at  the  current  “front”.  It 
follows  from  the  earlier  discussion  regarding  thresholds  that  at  the  current  “front”,  there 
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Figure  2-7.  Simulation  of  data  in  shown  in  Figure  1-6.  Top:  The  chemoattractant  profiles 
for  cells  a,  b,  and  c,  respectively.  Bottom:  The  dotted  lines  show  the 
quasisteady  state  profiles  attained  in  response  to  pipette  1.  The  full  lines  show 
the  new  quasisteady  state  profiles  attained  in  response  to  pipettes  1 and  2. 

Cell  a responds  to  the  new  chemoattractant  source  by  translocating  the 
phosphoinositide  peak  in  a wave-like  motion.  Cells  b and  c are  essentially 
unaffected  by  the  new  chemoattractant  source. 


37 


is  a tendency  for  the  threshold  to  increase  due  to  elevated  inhibitor  concentrations.  If  the 
perturbation  is  large  enough  at  the  current  “front”,  this  threshold  is  mitigated  and  a new 
peak  forms.  However,  at  the  “previous”  front,  the  tendency  of  the  threshold  to  increase 
due  to  elevated  inhibitor  concentrations  is  further  exacerbated  by  the  lower  active 
receptor  concentration.  The  thresholds  at  the  “previous”  front  become  so  large  that 
despite  the  large  concentrations  of  membrane  phosphoinositides,  they  fall  short  of  the 
threshold,  and  the  pre-existing  peak  collapses. 

Along  the  same  lines,  the  experiment  shown  in  Figure  1-7  can  be  simulated.  The 
model  cells  are  exposed  to  pipette  1 and  are  allowed  to  form  quasi-steady 
phosphoinositide  peaks  consistent  with  the  external  gradients  imposed  on  them  (dotted 
lines  in  Figure  2-7).  According  to  the  model  the  closer  the  cell  is  to  the  pipette,  the  more 
intensely  polarized  it  becomes  (dotted  lines,  Fig  2-7).  Subsequently,  all  three  cells  are 
exposed  to  an  additional  pipette  2 and  allowed  to  reach  a new  quasisteady  state  (full  lines 
in  Figure  2-7).  The  results  of  the  simulation  are  consistent  with  the  experimental 
observations 

• The  leading  edge  of  cell  a senses  a shallow  gradient  from  pipette  2.  However,  since 
pipette  1 exerts  a weak  influence  on  this  cell,  the  pre-existing  peak  is  easily 
dislocated  and  moves  like  a traveling  wave  towards  pipette  2 (Figure  2-7, Cell  a). 

• The  trailing  edge  of  cell  b experiences  the  strong  influence  of  pipette  2.  But,  the 
pre-existing  influence  of  pipette  1 is  so  strong  that  the  cell  does  not  entertain  the 
signal  from  pipette  2 (Figure  2-7,  Cell  b).  To  understand  this,  it  suffices  to  note  that 
the  stronger  the  pre-existing  polarity,  the  higher  the  concentration  of  the  inhibitory 
pool  of  I . Hence,  the  strength  of  the  signal  required  to  provoke  a new  peak 
increases  with  the  intensity  of  the  pre-existing  polarity. 

• The  leading  edge  of  cell  c senses  only  a weak  gradient  from  pipette  2.  In  addition, 
since  pipette  1 exerts  a strong  influence  on  cell  c,  the  phosphoinositide  peak 
remains  unaffected  (Figure  2-7c,  Cell  c). 
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Slow  Responses 

To  study  the  effect  of  adaptation  on  the  phosphoinositide  dynamics,  the  reduced 
model  (2.7-2.14)  should  be  considered. 

Adaptation 

To  simulate  the  effect  of  adaptation,  imagine  that  a resting  cell  is  suddenly  exposed 
to  a uniform  chemoattractant  concentration,  A+  . Since  the  environment  is  uniform, 
diffusion  plays  no  role,  and  the  evolution  of  the  system  is  given  by  ODEs  obtained  when 
8p  -8Pi  = 8j  = 0 in  equations  (2.7-2.14). 
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Figure  2-8.  Adaptation  in  response  to  a uniform  increase  in  the  chemoattractant 

concentration,  (a)  Evolution  of  the  adaptation  subsystem.  The  fast  variable, 
ju , rapidly  increases  from  zero  to  A7(l  + A+) , but  returns  to  zero  as  the  slow 
variable,  3 , gradually  increases  from  zero  to  /T/(l  + A+) . (b)  Evolution  of  the 
polarization  subsystem.  All  three  variables  change  rapidly  on  a fast  time  scale 
and  return  to  their  resting  levels  on  a slow  time  scale. 

Figure  2-8  shows  the  evolution  of  the  model  variables  for  A+  = 0.5 . Consistent 

with  the  experiments,  fu  and  n increase  on  a time  scale  of  1 second.  On  a slow  time 

scale  of  1000  seconds,  3 approaches  p(A+)  and  /j  approaches  zero,  so  that  n , ns  and 


i return  to  their  resting  values. 
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This  simulation  can  be  understood  by  appealing  to  Figure  2-8b.  In  the  absence  of 
diffusion,  the  non-uniform  steady  state  does  not  exist  and  the  uniform  steady  state  is 
stable.  The  system  therefore  evolves  along  the  locus  of  uniform  steady  states  (dotted 
curve  in  Figure  2-8b).  Initially,  n jumps  rapidly  from  its  basal  value  to  its  maximum 
value.  This  rapid  transient  is  followed  by  a slow  decline  along  the  path  of  uniform  steady 
states. 

Spontaneous  polarization 

The  above  simulation  of  adaptation  appears  to  contradict  the  experimental  data. 
For,  the  simulations  imply  that  the  cells  remain  in  a uniform  steady  state  throughout  the 
transient,  but  the  experiments  show  that  the  cells  frequently  polarize  before  they  adapt. 
This  contradiction  can  be  resolved  by  observing  that  single  cells  receive  noisy  signals. 
Indeed,  the  random  location  of  the  polarity  in  spontaneous  polarization  implies  that  some 
variable  that  is  upstream  of  the  phosphoinositides  undergoes  stochastic  fluctuations.  The 
most  upstream  source  of  stochastic  fluctuations  is  the  concentration  of  the 
chemoattractant.  We  show  below  that  fluctuations  of  the  chemoattractant  concentration 
are  sufficient  to  trigger  spontaneous  polarization  of  the  cell. 

To  model  the  fluctuations,  we  assume  that  the  uniformly  distributed 
chemoattractant  around  the  cell  is  in  thermodynamic  equilibrium.  Then  it  follows  from 
the  theory  of  fluctuations  [64]  that  the  chemoattractant  concentration  has  a Gaussian 
distribution  with  mean,  say  lm , and  variance  (, KkT/V)mlm , where  k is  the  isothermal 
compressibility  at  temperature  T , k is  Boltzmann’s  constant,  and  V is  the  characteristic 
volume  of  chemoattractant  sampled  by  the  receptors.  A reasonable  approximation  for  the 


40 


extracellular  sample  volume,  V , is  (Z),/r  )3/2  where  D,  is  the  diffusion  coefficient  of  the 


chemoattractant,  and  k is  the  dissociation  rate  constant  for  receptor-ligand  binding  [65], 
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Figure  2-9.  Spontaneous  polarization  in  response  to  a uniform  but  noisy  chemoattractant 
profile.  Parameter  values  used:  « = 10,  T = 298  K,  /c  = 45xlO-11  1/Pa, 

D,  = 10“10  m 2 Is.  (a)  Snapshot  of  the  noisy  chemoattractant  profile  generated 
by  the  stochastic  model,  (b)  Uniform  increase  and  return  to  rest  state  without 
any  spontaneous  polarization.  Observed  for  = 0.2.  (c)  Uniform  increase 
followed  by  the  formation  and  disruption  of  a phosphoinositide  peak. 
Observed  for  A+  = 0.45 . (d)  Uniform  increase  and  decrease  of 
phosphoinositides,  followed  by  peak  formation  and  disruption.  Observed  for 

r =1.5. 


If  we  partition  the  sample  volume  ( V ) into  n equal  sections,  the  evolution  of  the 
ligand  concentration  in  each  section  is  given  by  the  stochastic  equation 
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\rlcT 

Il(t)  = lm(l+l—zl),i  = l,2...n  (2.22) 

where  Vt  is  the  volume  of  each  section,  V In  and  z(.  is  a Gaussian  random  number 
generator  with  zero  mean  and  unit  variance.  Introducing  the  dimensionless  quantities 
X = l/(k+/k~ ) and  fa  = (/ ckT/Vj )1/2 , we  may  rewrite  equation  (2.22)  as 

Mr)  = Xm(\  + fazl),i  = \,2...n  (2.23) 

To  simulate  the  experiments  showing  spontaneous  polarization,  it  is  assumed  that  at 
r < 0 , the  cell  is  exposed  to  a small  and  uniform  chemoattractant  concentration,  Xm  = X" . 

At  t > 0 , the  cell  is  immersed  in  a mean  concentration,  X+  > X" . Sampling  the  random 
number  generator  (z.)  for  each  section  independently,  at  every  dr  = 0.01 , equation 

(2.23)  yields  a noisy  active  ligand  distribution.  A typical  snapshot  of  this  distribution  is 
shown  in  Figure  2-9a.  The  response  to  the  noisy  ligand  concentration  is  then  simulated  by 
introducing  X(t,£)  in  the  model  equations  (2.7-2.14). 

Depending  on  the  magnitude  of  X+ , three  types  of  responses  are  obtained  (Figure 
2-9b,  c,  d).  These  responses  are  consistent  with  the  experimental  observations  of  [66],  To 
understand  these  responses,  observe  that  the  noise  in  the  ligand  concentration  induces 
minute  gradients.  Beyond  a critical  X*ean , this  non-uniform  environment,  makes  the  rest 

state  unstable  in  the  Turing  sense,  i.e.,  stable  to  homogenous  perturbations  and  unstable 

to  the  smallest  non-homogeneity  [67].  Hence,  as  the  system  adapts,  it  evolves  along  the 

path  of  nonuniform  steady  states  on  the  bifurcation  diagram  (Figure  2-5b).  Thus 

1.  If  X+  is  small,  n increases  uniformly  and  then  returns  to  its  rest  state  without  any 
spontaneous  polarization  during  the  transient  (Figure  2-9b). 
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2.  At  intermediate  values  of  A+ , n increases  uniformly  on  a fast  time  scale,  but  then 
breaks  out  into  a nonuniform  distribution.  As  time  evolves,  the  width  of  the  peak 
progressively  decreases  until  the  peak  disappears  completely  and  the  system 
returns  to  its  rest  state  (Figure  2-9c). 

3.  If  A+  is  large,  n increases  uniformly  on  a fast  time  scale,  and  then  decreases 
uniformly  on  a slow  time  scale  before  erupting  into  a nonuniform  distribution. 

The  polarization  gradually  decreases  and  the  system  returns  to  its  rest  state 
(Figure  2-9d). 

Here,  adaptation  of  variable  fi  brings  the  system  back  to  its  rest  state  and  causes 
disruption  of  the  peak. 

These  results  are  consistent  with  the  sensory  adaptation  observed  in  human 
neutrophils  [29],  The  time  taken  for  the  rapid  uniform  increase  is  relatively  independent 
of  the  chemoattractant  concentration.  On  the  other  hand,  the  time  taken  for  the 
spontaneous  polarization  to  occur  increases  with  the  mean  chemoattractant  concentration. 


CHAPTER  3 

MATHEMATICAL  ANALYSIS 

Since  the  molecular  basis  of  gradient  sensing  is  not  completely  known,  it  seems 
premature  to  expect  quantitative  answers  from  the  mathematical  model.  However,  one 
can  still  appeal  to  the  mathematical  model  for  answers  concerning  qualitative  properties. 
For  example,  it  is  quite  relevant  to  ask:  What  are  the  factors  that  determine  the  geometry 
of  the  phosphoinositide  localization?  What  are  the  conditions  that  allow  for  the  formation 
of  the  phosphoinositide  localization  from  a homogenous  resting  state?  Such  questions  can 
be  answered  by  a mathematical  analysis  of  the  model. 

The  polarization  subsystem  of  the  model  is  a boundary  value  problem  with 
periodic  boundary  conditions  and  the  phosphoinositide  localization  is  a non-uniform 
steady  state  solution  of  this  problem.  Since  the  diffusivity  of  membrane  phospholipids  is 
three  orders  of  magnitude  smaller  than  the  diffusivity  of  the  cytosolic  pool  of  inositol 
phosphates,  the  boundary  value  problem  contains  a small  parameter,  namely  the  ratio  of 
the  two  diffusivities.  Approximations  to  the  non-uniform  steady  states  can  therefore,  be 
obtained  by  appealing  to  a technique  called  singular  perturbation  theory.  This  approach, 
originally  developed  by  [68],  has  been  successful  for  analyzing  activator-inhibitor  models 
in  which  the  inhibitor  diffuses  much  more  rapidly  than  the  activator.  The  first  half  of  this 
chapter  describes  this  approach  for  the  analytical  construction  of  the  phosphoinositide 
steady  state.  It  will  allow  us  to  obtain  approximate  geometrical  properties  of  the 
phosphoinositide  steady  state  such  as  width  and  height.  We  shall  also  use  this  study  to 
understand  how  variation  in  phosphoinositide  synthesis  can  lead  to  change  in  peak  width. 
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The  second  half  of  this  chapter  deals  with  a stability  analysis  for  the  homogenous 
phosphoinositide  resting  state.  This  analysis  provides  insight  into  the  range  of 
phosphoinositide  synthesis  rates  that  cause  the  homogenous  resting  state,  to  lose  stability 
and  spontaneously  polarize  into  a phosphoinositide  localization. 

Simplified  Polarization  Subsystem 

It  has  been  experimentally  shown  that  the  addition  of  extracellular  membrane 
permeant  PIP3  to  resting  cells  can  initiate  spontaneous  polarization  of  intracellular 
phosphoinositides,  morphological  polarity  and  motility,  in  the  absence  of  any  added 
chemoattractant.  This  suggests  that  the  polarization  subsystem  can  be  stimulated 
independent  of  receptor  activation  and  the  adaptation  subsystem. 

This  experiment  can  be  modeled  by  perturbing  the  intracellular  concentration  of 
phosphoinositides  and  studying  the  response  of  the  polarization  subsystem.  The  model 
can  be  further  simplified  by  assuming  that  the  sum  of  phosphoinositide  concentrations  in 
the  membrane  and  endoplasmic  reticulum  is  a constant  (pt  = p + ps  = constant ),  over  the 

time  scale  of  the  experiment.  This  assumption  reduces  the  polarization  subsystem  to  the 
following  equations, 


~ = kfp2(p  - p)-k  pi  + c -k  p-\ — 
dt  ft?  yj  ry  P pf  R2  dez 


di 

~dt 


~s{kfp2  (p,  -p)~  krpi)  + c,~  k,i  + -i 


D,  d2i 


R 2 dd1 


with  periodic  boundary  conditions, 


x = p,i 


(3.1) 

(3.2) 


(3.3) 


and  initial  conditions 
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p{e,Q)  = p-{6),  i(0,O)  = r(0),O<0<2;r  (3.4) 

where  p~{6)  and  i~  (9)  are  chosen  to  reflect  the  perturbation  imposed  in  the  experiment. 
Using  the  dimensionless  variables 


P 1 c 

n=~,i  = — 


e 


,T  = - 


p,  spt  2x3.1416  \l{krsp') 


and  dimensionless  parameters 


8 =i^- 

u p Kv, 


D,/C 2 


t krsp, 


cjp, 

cjp. 


kB 

ic  = p— 

p “ K*Pt 


i Ksp, 


K.  = 


Ksp, 


equations  (3.1)-(3.4)  can  be  rewritten  as 


dn  2 dn 

— = Krn  {\-n)-m  + u/„-K„n  + 5 — - 
Sr  f r r p d£2 


A r, 

Sr  f 

with  dimensionless  boundary  conditions 


di  2n  x o d2i 

= Kfn  (1  -n)-ni  + if/i-Kji  + oi- 


S£2 


x(0,r)  = *(U),^  = ^U>0 
S^  3£ 

X = 7t,l 


and  dimensionless  initial  conditions. 


(3.5) 

(3.6) 


(3.7) 


^(0,£)  = 7T~,i(0,^)  = T,0  < £ < 1 (3.8) 

fc f — 1.5  has  been  used  for  the  following  simulation.  All  other  parameter  values  are  the 

same  as  in  Table  2-1.  Before  the  phosphoinositide  perturbation  is  applied  (r  < 0),  P and 
I are  in  a homogeneous  steady  state  {n , f ),  determined  by  the  equations 


Kf  1 - n)  - ni  +y/p-  Kpn  = 0 

*7;r2(l -n)  -m  +ipl  - k,J  = 0 
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Figure  3-1.  Steady  states  of  the  reduced  model.  Dashed  lines  denote  initial  conditions; 
Solid  lines  denote  steady  state  response,  (a)  Phosphoinositides  %)  (b) 
Inositol  phosphates  i(£) . 

At  t = 0,  the  homogenous  steady  state  of  the  phosphoinositides  (n)  is  subjected  to 
a non-uniform  perturbation  ( n~  = n + A;r(£) ) such  that  concentration  is  highest  at  the 
leading  edge  (£  = 1/2 ) and  lowest  at  the  trailing  edge  (£  = 0,1 ).  A pronounced 


phosphoinositide  peak  develops  at  the  leading  edge  (£  = 1/2 ) (Figure  3-la).  Compared  to 
the  polarized  distribution  of  membrane  phosphoinositides,  the  concentration  profile  of 
inositol  phosphates  is  essentially  flat  (Figure  3- lb).  All  the  analysis  performed  in  the 
following  sections  will  be  based  on  this  simplified  model. 

Evolution  of  the  Phosphoinositide  Localization 
The  evolution  of  the  phosphoinositide  localization  is  a two-fold  process.  Firstly,  if 
a certain  threshold  is  crossed,  membrane  phosphoinositides  develop  large  spatial 
gradients  from  their  homogenous  initial  state.  This  evolution  reaches  a state  of  quasi 
equilibrium  rapidly.  In  the  second  phase,  the  large  gradients  that  form  move  outward  on  a 
slow  time  scale,  until  a steady  state  is  reached. 
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Figure  3-2.  Nullcline  dn / dx  — f(n,i)  — 0 . All  points  above  the  curve  (dn  / dx  <0)  move 

towards  lower  n and  all  points  below  the  curve  (dn ! dx  >0)  move  towards 
higher  n . 

First  Phase 

Since  the  spatial  gradients  in  P are  initially  small,  diffusion  can  be  neglected  and 
the  initial  evolution  of  P can  be  approximated  by 


dn 

dx 


Krn 2 (1  - *)  - m +^P~  Kpn  = —f(n,  i ) 
0)  = n~  = n + A?r(£) 


(3.9) 


where  n(<^,0)  denotes  the  initial  distribution  of  P . During  this  period,  it  is  assumed  that 

the  concentration  of  i is  held  fixed  at  its  initial  value  f = f . If  n at  the  leading  edge  of 
the  cell  crosses  the  threshold  for  amplification,  i.e.  if  the  perturbation  is  in  the  region 
where  ^->  0,  the  initial  point  immediately  increases  towards  higher  n (see  Figure  3-2). 
In  contrast,  if  the  trailing  edge  of  the  cell  does  not  feel  the  perturbation,^  will  lie  below 
the  threshold  for  amplification  (the  region  where  < 0 , see  Figure  3-2),  forcing  n to 
either  decrease  or  stay  at  n = n~  . 
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Since,  the  phosphoinositides  in  the  membrane  either  reach  a high  steady  state  or  a 
low  steady  state,  the  first  phase  of  evolution  results  in  a state  of  quasi-equilibria  - two 
smooth,  distinct  subregions  separated  by  sharp  transition  layers. 

Second  Phase 

At  the  end  of  the  first  phase,  each  of  the  subregions  reaches  a steady  state. 
However,  the  sharp  transition  layers  between  these  subregions  are  in  a state  of  imbalance, 
since  there  is  more  synthesis  of  phosphoinositides  than  removal.  The  imbalance  in 
phosphoinositides  causes  the  walls  of  the  transition  layers  to  move  outward  like  a 
traveling  wave,  until  there  is  no  more  accumulation  of  phosphoinositides  within  the 
transition  layer.  This  eventually  results  in  the  formation  of  the  steady  state, 
phosphoinsoitide  localization. 


with  periodic  boundary  conditions  //(0)  = //(l)  and (0)  =-J(l) , where  ii  = n,i. 
Here  f(n,i ) denotes  the  rate  of  removal  of  membrane  phosphoinositides  ( P ) and 
g{n,i)  the  rate  of  removal  of  cytosolic  inositol  phosphates  (I). 


Construction  of  the  Non-homogenous  Steady  State 

The  non-homogenous  steady  state  of  the  model  satisfies 


(3.10) 


(3.11) 


As  discussed  earlier,  it  is  important  to  note  that  a stable  non-homogenous  steady 
state  in  n can  form  only  when  Sp  <K  8t . Problems  of  this  type  where  a small  parameter 

8p  multiplies  the  highest  derivative  usually  give  rise  to  singular  solutions,  that  undergo 
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rapid  transitions  in  £ [69],  In  these  cases,  it  becomes  difficult  to  find  a complete 
analytical  solution  that  is  valid  for  all  0 < £ < 1 . The  standard  way  to  tackle  this  difficulty 
is  to  seek  a solution  which  ignores  diffusion  in  n everywhere  except  for  the  region  where 
n changes  sharply.  In  other  words,  it  is  assumed  that  the  solution  comprises  of  two 
regions:  (1)  An  ‘outer’  solution  where  n changes  slowly  with  £ and  diffusion  can  be 
neglected  (2)  An  ‘inner’  solution  where  n changes  fast  with  £ and  diffusion  needs  to  be 
taken  into  account.  The  outer  solution  will  yield  smooth  subregions  that  are  separated  by 
discontinuities,  which  arise  because  diffusion  of  P has  been  ignored.  The  inner  solution 
corresponds  to  the  sharp  interface  that  should  exist  at  these  points  of  discontinuities.  On 
patching  up  these  two  solutions,  an  approximate  steady  state  peak  can  be  constructed  [67, 
70,71], 

Outer  Solution 

First,  assume  that  the  diffusion  of  n is  small  everywhere  (i.e.  Sp  ^ * 0 ).  Setting 
S,  - 1 and  using  regular  perturbation  expansions  in  small  parameter  Sp , 

n = 7r(0)  + 8pnm..  and  i = i(0)  + Spim  ..,  the  zeroth  order  equations  ( Sp  ~ 0 ) for  (3.10)  and 

(3.11)  can  be  written  as 


The  function  f(n,  i)  = -Kfn2  (1  - n)  + m - yp  + Kpn  = 0 is  an  S-shaped  curve  that 
has  a maximum  in  i at  ( n0 , i0  ) and  a minimum  at  ( n0 , i0 ) (Figure  3-2,  3-4a). 


f(n,  i)  = -icf7r2(l  - n)  + ni  — \yp  + Kpn  = 0 


(3.12) 


(3.13) 
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Figure  3-4.  Phase  portraits  for  each  branch  of  / {n,  i)  = 0 (a)  Nullclines  / {n,  i)  = 0 and 
g(n,i)  = 0 . For,  ^ = 1.5,  g(n,i)  = 0 intersects  the  lower  branch  of 
f{n,i)  = 0 at  (nh,ih)  (b)  Phase  portrait  for  lower  branch  I (c)  Phase  portrait 
for  middle  branch  II  (d)  Phase  portrait  for  upper  branch  III 

The  nullcine  f{n,i)-  0 has  three  separate  roots 


n(i)  m < 


n,  (0 

n (t) 

n v 7 

n 0) 

in  v 7 


If  n <7t0  - Lower  Branch  I 
If  n Q < n < nQ  - Middle  Branch  II 
If  n > 7r0  - Upper  Branch  III 


(3.14) 


Using  the  roots  of  f(n,i)  = 0 , equation  (3.13)  can  be  reduced  to 
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d2i 


= g(Xj(0,i)  wherey  =1,  II,  III 


(3.15) 


with  periodic  boundary  conditions  z(0)  = i(l)  and  ^(0)  = -^(1) . Since  the  diffusion  of  I 
is  fast,  only  continuous  solutions  for  *(£)  are  expected.  In  order  to  solve  for  the  outer 
solution,  z(£),  the  equation  (3.15)  can  be  rewritten  in  the  form  of  a conservative  system, 

^-  = ±J2(E -</,,)  (3.16) 


where  $ denotes  the  potential  of  the  system  and  E(i,i  'f)  is  the  Hamiltonian. 

Equation  (3.16)  can  be  used  to  draw  the  phase  plane  (Figure  3-4),  for  each  branch 
of  n{i) . Any  orbit  that  starts  from  a point  on  the  phase  plane  (at  £ = 0 ) and  returns  to  the 
very  same  point  within  £ = 1 would  be  a valid  solution  of  (3.15)  and  would  satisfy  the 
periodic  boundary  conditions. 

Branch  II  can  yield  valid  solutions  to  the  problem  (3.15)  (Figure  3-4c).  However, 
solutions  arising  from  this  branch  turn  out  to  be  unstable,  and  therefore  have  to  be 
dismissed  [67].  From  inspecting  the  orbits  of  branch  I and  branch  III  (Figure  3-4  b,  d), 
one  can  conclude  that  none  of  these  orbits  satisfy  periodic  boundary  conditions  - since 
none  of  the  orbits  close. 

Closer  examination  of  these  phase  portraits  reveals  another  possibility:  A solution 
with  the  required  period  (see  figure  3-5,  orbit  ABCB’A)  can  be  constructed  by  patching 
together  orbits  of  branch  III  with  the  orbits  of  branch  I.  Since,  every  point  on  this  unique 
orbit  can  be  chosen  as  a starting  point  (£  = 0) , infinite  such  solutions  are  possible.  These 
solutions  are  just  lateral  translations  of  one  another  in  space. 
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For  specificity,  let  us  restrict  ourselves  to  constructing  a solution  that  has  one 
maxima  (n=l)  at  % = 1/2.  Solutions  with  more  than  one  maxima  (n>l)  can  also  be 
constructed  (see  n=2.  Figure  3-5).  However,  these  solutions  turn  out  to  be  unstable  [67]. 


Figure  3-5.  Patched  phase  plane  for  the  outer  solution i{%) . The  outer  orbit  ABCB’A  is  a 
solution  with  a single  maximum  in  i(£)  (n=l).  The  inner  orbit  AC  is  a 
solution  that  has  two  maxima  in  i( £)  (n=2). 

In  order  to  construct  the  (n=l)  solution  (ABCB’A),  assume  an  arbitary  point, 
t0<i‘  <i0,  for  patching  up  branch  I and  branch  III.  Consider  an  orbit  that  starts  at  a point 
A(a,0)  (with  a >i0)  on  the  branch  I phase  plane  (see  Figure  3-5)  and  follow  the  orbit  to 
a point  B (i*,;‘ ) . Transfer  this  point  to  the  branch  III  phase  plane,  and  follow  the 
corresponding  orbit.  If  this  orbit  reaches  point  C(/3, 0)  in  £ = 1/2 , the  orbit  would  satisfy 
(3.15)  and  its  boundary  conditions.  Note  that  only  a specific  choice  of  a will  pick  up  the 
desired  solution.  This  particular  a can  be  found  by  imposing  (i)  continuity  of  the 
patched  solution  at  i = t and  the  (ii)  size  of  the  system  ( A£  = 1) . 
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In  order  to  patch  up  the  orbits,  let  us  define  the  Hamiltonian  of  a particular  orbit 

= (3.17) 

l 

where  the  potential  <j)(‘Jn)  =-^g{n,  m{i),i)di . 

The  solution  to  (3.15)  must  satisfy  the  following  conditions. 

1 . Continuity  at  i = C : In  order  to  ensure  the  continuous  transition  of  the  branch  I 
solution  to  branch  III  and  vice  versa  (at  B and  5',  Figure  3-5a)  the  following 
conditions  have  to  hold, 


1(0 

, \W) 

~li  ~l( 

(3.18) 

!(/) 

(3.19) 

■l  ~ 

=l*l 

l (,) 

BJ? 

=i<m)=r 

BJ? 

(3.20) 

Using  (3.17),  (3.20)  and  the  fact  that  the  Hamiltonian  of  a particular  orbit  is 


constant,  the  following  condition  should  hold  for  branch  I 

EV\i,i()  = E«\a,0)  = E^  (3.21) 

=>#a)=^(i'f|(/>)2+^/)  (3-22) 

2 BJ* 

Similarly,  for  branch  III 

£w(Mf)  = E(cm)(/3, 0)  = E™  (3.23) 

<3.24) 

2 B 

Using  (3.18),  (3.19),  (3.22)  and  (3.24)  the  following  condition  can  be  written 

J!  g(*  (0 ,i)di  = Jf  g{nm  ( i),i)di  (3.25) 


54 


2.  Periodic  Boundary  Conditions  at  ^ = 0,1:  The  time  taken  for  moving  through  the 
orbit  ABCB’A  should  be  equal  to  the  width  of  the  system  ( = 1).  This  can  be 

used  to  derive  a time-map  equation.  Rearranging  (3.17)  for  branches  I and  III,  and 
using  (3.21),  (3.23)  we  have 


Here,  and  denote  the  £ values  at  points  B and  B’.  Adding  these  four 
integrals,  the  following  time-map  condition  can  be  written 


For  a given  i - 1 , a and  (3  can  be  found  by  solving  conditions  (3.25)  and  (3.28) 
simultaneously.  Knowing  or  and  /?,  the  unknowns can  be  found  using 

equations  (3.22)  and  (3.27).  The  unknown  patching  point,  i = i can  be  determined 
from  a continuity  condition  for  the  inner  and  outer  solution  of  the  problem.  This 
will  be  discussed  in  the  following  section.  Assuming  i is  known,  the  following 
equations  can  be  numerically  solved  for  the  outer  solution  z(£), 


(3.26) 


Using  (3.26)  the  following  equations  can  be  written 


(3.27) 


(3.28) 


(3.29) 


f(^j,i(l))  = 0 =>  it  — nj(iU)) 


i =uu 


(3.30) 
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using  (i)  Periodic  boundary  conditions  on  *(/)( £)  at  £ = 0,1,  i{m)(\/2)  = J3  and 

i4u,,)(\I2)  = 0 (ii)  Continuity  conditions  (3.18),  (3.19)  and  (3.20).  Although  i is 

continuous  at  £ = %BB.,  the  corresponding  solution  for  ;r(£)  (from  (3.30))  has  a 

discontinuous  transition  from  branch  I to  branch  III  or  vice  versa  (see  Figure  3-6), 
i.e., 

*(#-,)  = *(#+,)**,(*•) 

<3J,) 

Here,  £ refers  to  a small  displacement  on  either  side  of  £ = £ . 

BJ¥ 

Hence,  the  outer  solution,  ;r(£)  (Figure  3-6),  consists  of  two  distinct  subregions, 
a low  steady  state  *>(*(£))  and  a high  steady  state  ^ (/(£)) . As  expected,  each  of  these 
subregions  is  separated  by  discontinuities  at  points  of  patching  £ = £B  s, . The 

concentration  of  n in  each  of  these  subregions  alters  itself  according  to  the  local 
concentration  of  i , such  that  there  is  no  net  formation  of  membrane  phosphoinositdes 

(/0m)  = 0). 

The  above  calculations  also  allow  us  to  evaluate  approximate  geometrical 
properties  of  the  phosphoinositide  localization,  such  as  width  of  the  localization  A -E 
and  height  of  the  localization,  *,„(/(£*  *.))-*,  (*(&.„.))  - 
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Figure  3-6.  Outer  solution  n{^) . Discontinuities  exist  at  £B  and  £B, . Solution  between 
and  £B,  corresponds  to  nlu  (i(g)) . Solution  to  the  left  of  and  to  the 
right  of  £B,  , corresponds  to  n,  («(£)) . 

Therefore,  at  the  end  of  the  first  stage  of  evolution,  the  phosphoinositides  reach  an 
upper  steady  state  nm(i(^)) , at  the  site  of  the  perturbation  and  a lower  steady  state 

nm  (i(£)) , at  all  other  regions  of  the  membrane.  The  interface  that  should  exist  at  the 

point  of  discontinuity  corresponds  to  the  inner  solution  to  the  problem. 

Inner  Solution 

In  the  previous  section,  Sp  was  assumed  negligible.  This  assumption  yielded  a 
discontinuity  in  n , at  £ = <^B  B, . However,  the  true  solution  to  the  problem  (for  5 * 0) 
actually  has  a sharp  interface  in  this  region  (see  Figure  3- lb).  This  discrepancy  occurs 
because  8p  ^ has  been  neglected  in  the  region  of  sharp  change.  A regular  perturbation 

analysis  does  not  pick  up  the  singular  change  near  £ = £BjB.  . In  fact,  it  shrinks  the  region 
to  a point. 

In  order  to  handle  this  difficulty  another  solution  needs  to  be  found  near  £ = £B  fl, . 
This  can  be  done  by  stretching  the  £ coordinate  in  the  region  and  imposing  boundary 
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conditions  that  the  solution  should  approach  the  outer  solution  as  we  move  on  either  side 
of  E,  — £B'B. . On  rescaling  the  region  near  £ = E,B  B.  as 


y=- 


S-Sb, 


& 


equations  (3.10)  and  (3.1 1)  can  be  reduced  to 

d2n' 


dy 2 


= /(*',*') 


2.i 


2 =sPg(*,j) 


(3.32) 


(3.33) 

(3.34) 


Here  the  superscript  ‘ / ’ denotes  the  inner  solution  to  the  problem.  Using  regular 
perturbation  expansions  in  small  parameter  8p  for  n‘  =n[+8pn[..  and  i‘  =i‘0  + 8pi[..,  the 
zeroth  order  equations  for  (3.33)  and  (3.34)  can  be  written  as 


d2n‘ 

dy 2 


= /(*',*') 


d2i‘ 

dy2 


-o 


(3.35) 

(3.36) 


Note  here,  that  the  small  term  on  this  length  scale  is  8p  g(;r',f)and  not  the  second 
derivative  term.  The  denominator  of  (3.32),  ^8~p  , is  a small  quantity.  Therefore,  a large 
change  in  y is  required  to  shift  £ on  either  side  of  ^B  B, . Moving  slightly  to  the  left,  say 
£ = # - . translates  to  y -»  -go  . Similarly,  E - E + translates  to  y ->  oo . 

Now,  consider  the  region  around  E,  = £ . Continuity  of  the  zeroth  order  inner  and 
outer  solution  ensures  that 


7r'(-°o)  = 7rf  0*) 


(3.37) 
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7r‘(+co)  = 7rm  (O 


(3.38) 


l'(-0o)  = f(+O0)=/‘ 


(3.39) 


Using  (3.36)  and  boundary  condition  (3.39),  we  can  get  i'(y)  = i . 

The  unknown  parameter  i can  be  found  using  the  continuity  conditions  (3.37) 
and  (3.38).  Multiplying  (3.35)  by  integrating  from  -oo  to  +ooand  using 

^■(i°°)  — > 0 , we  get 


Equation  (3.40)  has  the  following  physical  interpretation.  During  the  second  stage 
of  evolution,  the  sharp  interface  moves  outward  like  a traveling  wave  due  to  an 
accumulation  of  phosphoinositides.  As  the  wave  propagates  outward,  the  local 
concentration  of  the  inositol  phosphates  («’ ) keeps  changing.  When  the 
phosphoinositides  within  the  interface  are  produced  at  the  same  rate  as  they  are  removed 
(i.e.,  when  (3.40)  is  satisfied),  the  wave  stops  propagating  and  forms  a stationary  front. 
This  balance  in  phosphoinositide  synthesis  and  removal,  determines  the  width  of  the 
phosphoinositde  localization. 

In  order  to  solve  the  inner  solution  to  the  problem,  (3.35)  can  be  rewritten  in  the 

form 


+oo 


(3.40) 


-oo 


dn 


(3.41) 


where  <f>x  = -^f{n‘,i)dn  is  the  potential  of  the  system  and  E{n\n'‘y)  is  the 


Hamiltonian.  Using  (3.41),  we  can  draw  the  phase  plane  of  (3.35),  for  / = /* . Any  orbit  on 
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this  phase  plane  that  satisfies  (3.40)  and  the  boundary  conditions  (3.37),  (3.38)  is  a valid 
inner  solution  to  the  problem.  Equation  (3.40)  implies  that  d>  . = 6 . . This  in  turn 

implies  that  Efn^i), 0)  = E(nm{T), 0) . The  only  orbit  that  can  satisfy  all  these  above 
conditions  is  a heteroclinic  orbit,  7r'~(y) , (AB,  Figure  3-7)  that  connects  the  two  saddle 
points,  n,{i)  and  nm(i).  Similarly,  the  inner  solution  at  %B,  corresponds  to  the 
heteroclinic  orbit,  n,+{y) , (BA,  Figure  3-7)  that  connects  n a (/*)  to  n (f). 


Figure  3-7.  Heteroclinic  orbit  (AB)  representing  inner  solution  n\y). 

Zeroth  Order  Approximation  to  Complete  Solution 

The  zeroth  order  solution  (see  Figure  3-8)  to  the  polarization  subsystem  (3.1),  (3.2) 
can  be  found  by  patching  up  the  inner  and  outer  solutions.  Using  the  fact  that  the  inner 
solution  asymptotically  approaches  constant  values,  or  nUJ{%BB.)  on  either  side 

of  the  sharp  interface,  the  zeroth  order  solution,  nc{^),  can  be  written  as 
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*c(f)  = 


<5P 

o*4*4B 

*w(4)+xi0  J-B)  *0,(4,) 

<SP 

4b - 4 - 1/2 

ylsp 

^*4*4V 

<5p 

4,’*4*  i 

(3.42) 


Figure  3-8.  Zeroth  order  solutions  (dashed  lines)  in  comparison  with  the  numerically 
calculated  complete  solution  (solid  lines),  (a)  Plot  of  . (b)  Plot  of  i( £) 
(corresponds  to  orbit  ABCB’A  in  Figure  3-5). 

Since  the  inner  solution  r'(£)  is  a constant, i(,\4BB.)  =r(///)(^BB.)  = «‘,  the  zeroth  order 
solution,  ic{%) , is  nothing  but  the  outer  solution  (Figure  3-8b) 


a#)= 


(4) 
r\t) 
pm 


4,*4*4g' 
4**4*  i 


(3.43) 
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Instability  of  the  Homogenous  Steady  State  with  Respect  to  Change  in  k} 

As  discussed  in  Chapter  2,  the  homogenous  resting  state  of  the  phosphoinositides 
evolves  into  a non-homogenous  localization  only  when  the  synthesis  rate  of  the 
phosphoinositides  (xy)  is  optimal.  This  property  of  gradient  sensing  can  be  analyzed  by 

performing  a Turing  stability  analysis  to  evaluate  the  stability  of  the  homogenous  steady 
state  to  non-homogenous  perturbations,  as  Kf  is  varied. 

Consider  the  uniform  resting  state  of  the  polarization  subsystem  i.e.  the 
homogenous  steady  state  obtained  when  no  stimulus  is  applied  to  the  model  cell,  n and 
i can  be  obtained  from  the  following  equations 


In  order  to  study  the  stability  of  the  resting  state,  we  may  impose  an  arbitrarily 
small  perturbation,  n +(£,r),f+(£,r)  to  fc  and  f , respectively.  If  the  resting  state  is 
stable, ;r  and  i will  return  to  their  homogenous  resting  state (n,i) . On  the  other  hand,  if 
the  resting  state  is  unstable,  u and  i will  leave  the  resting  state  and  evolve  into  another 
steady  state  in  long-time.  The  equations  that  describe  the  evolution  of  these  perturbations 
can  be  obtained  by  linearizing  equations  (3.5),  (3.6)  using  n = n + n+(%,T), 
i = i + i+(Z,T),  (3.44)  and  (3.45) 


f{n , f ) = -n)-  ni  + t//p-Kpn)  = 0 


(3.44) 


g{n , f)  = -ji)-  m +y/.  - Kft)  = 0 


(3.45) 


(3.46) 


(3.47) 
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where  f\,f\,g\,g\  denote  partial-derivatives  evaluated  at  the  resting  state 


Now,  the  solution  vector  u = (n+ ,i+)  for  (3 .46), (3, 47)  can  be  expanded  as  a 
Fourier  series  expansion 


u = 


n 


(3.48) 


where  <j>j{x)  denotes  the  jth  eigenfunction  of  the  eigenvalue  problem4 


d e 


= k2Mt) 


(3.49) 


with  periodic  boundary  conditions  ^(0)  = (1)  and  ^(0)  = ^(1) . Substituting  the 

expansion  (3.48)  into  (3.46)  and  (3.47),  and  equating  the  coefficients  for  each  ^(£) , we 
get  the  following  problem  for  s /(r) , 


ds,  _ 

— - = M,s,  where  M denotes 
dr 


J J 


-f, 


(3.50) 


S*  -kjd'-g., 

The  solution  to  the  linear  problem  (3.50)  should  be  of  the  form  c,v;1ev  +c2vj2eAjlt , 
where  c,,c2  are  constants,  , Xj2 denote  the  eigenvalues  of  M } and  xjX,xj2  denote  the 


corresponding  eigenvectors. 

Since,  each  matrix  M.  has  two  eigenvalues,  the  stability  of  the  homogenous 
steady  state  depends  on  2xj  eigenvalues,  ^ , A,2  etc.  If  any  of  these  eigenvalues 


4 The  eigenfunctions  of  the  operator  -d2  / d^2  have  the  form  <j>j  (£)  = a}e  tjt . Here,  kj  represents  non- 
negative eigenvalues  or  perturbation  wavenumbers,  2 jn  for  j=0,l,2... 
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have  positive  real  parts,  the  perturbation  grows  exponentially  and  the  homogenous  steady 
state  becomes  unstable.  On  the  other  hand  if  all  the  eigenvalues  have  negative  real  parts, 
the  perturbation  decays  and  the  homogenous  steady  state  is  stable. 

The  eigenvalues  for  each  M } satisfies  the  characteristic  equation 


V + W, + * ', + k)Sp  + k)S: ) + (k)5rs,  + k)S,  g ',+  k)S,f\ +/  V g',  - / = 0 

It  follows  that  the  stability  of  the  homogenous  steady  state  requires 


Trace(  Mj ) = -(/  \+ g \+ k)8p  + k)5i ) < 0 (3.51) 

Det(  Mj ) = (k)SpSt  + k)8pg  \ + k)8J \ + f\ g\- f\g\)>  0 (3.52) 
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Figure  3-9.  Plots  of  Det(Mo)  and  Trace(Mo)  as  parameter  k f is  varied. 

First,  consider  the  stability  of  the  homogenous  steady  state  with  respect  to 
spatially  homogenous  perturbations,  n+(j),i\r) . This  corresponds  to  setting k =0  and 
neglecting  diffusion  for  n , i . In  this  case  it  turns  out  that  as  Kf  is  varied,  the  conditions 


for  stability,  Trace(  M0  )<0  and  Det(  M0  )>0,  are  satisfied  for  all  tcf  (see  Figure  3-9).  It 
also  follows  that  there  is  only  one  unique  steady  state  {n,  f ) for  each  Kf . This  can  be 
shown  by  using  the  following  expressions 
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dft  _ (f\g\,-g\f'Kf) 
dKf  Det(M0 ) 

dr  (g'./y/uv,) 

dKf  Det(M0 ) 


(3.53) 


(3.54) 


Since,  the  denominator  of  (3.53), (3.54)  is  always  positive  for  all  Kf,  the  slope  of 
the  steady  state  curve  can  never  reach  infinity.  Hence,  as  Kf  is  varied,  only  one 


homogenous  steady  exists  and  this  steady  state  is  stable  with  respect  to  homogenous 
perturbations. 


axis  is  tangential  to  the  Det  (MJ)  curve.  Here,  kc  denotes  the  critical  wave 
number  at  which  instability  occurs,  (b)  Bifurcation  diagram  for  the 
homogenous  steady  state  as  Kf  is  varied.  Dashed  line  denotes  the  region  of 
Turing  instability;  Solid  line  denotes  stable  steady  states. 

Using  the  above  result,  the  Turing  stability  of  the  homogenous  steady  state  to 
non-homogenous  perturbations,  can  be  analyzed  using  the  conditions  (3.5 1),(3.52)  for  all 
k . First,  observe  that  condition  (3.51)  always  holds  as  Kf  is  varied.  This  follows  from 
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Trace(  Mf ) = Trace(M0)  - (k]Sp  + k2St)  < 0 (3.55) 

because  Trace(M0)  <0  and  (k2Sp+k2jSi)> 0. 

Hence,  instability  can  occur  only  via  condition  (3.52).  If  we  now  assume  that  k is 
continuous,  then  the  following  conditions  have  to  hold  at  the  limit  of  stability  (see  Figure 
3- 10a) 


Det(M;)  = 0 


dDet(Mj) 

~~dk 


= 0 


(3.56) 


Condition  (3.56)  implies  that  the  homogenous  steady  state  is  Turing  unstable  iff 

SpS\+S,f\  < -2  (Sps,  Det(M0))1/2  (3.57) 

This  instability  is  realized  with  respect  to  perturbations  whose  wave  number  (kcJ,  a 
discrete  wave  numbers  among  2}n  for  j=I,2...)  is  closest  to  kc  =(^,)‘1/2(Det(M0))1/4,  the 

critical  wave  number  at  which  condition  (3.57)  is  valid  (see  Figure  3-10a). 

Using  (3.57),  we  can  plot  the  bifurcation  diagram  for  the  homogenous  steady 
state,  as  Kf  is  varied  (Figure  3- 10b).  As  expected,  the  homogenous  steady  state  becomes 


Turing  unstable  only  at  intermediate  Kf.  Under  these  conditions,  even  the  smallest  non- 

homogenous  fluctuation  causes  the  homogenous  steady  state  to  lose  stability  and 
spontaneously  polarize  into  a non-homogenous  phosphoinositide  distribution. 


CHAPTER  4 

EXPERIMENTAL  ANALYSIS 

It  has  been  shown  in  Chapter  2,  that  the  model  simulations  are  consistent  with  the 
spatiotemporal  dynamics  observed  during  gradient  sensing.  The  specific  aim  of  this 
chapter  is  to  obtain  experimental  evidence  that  verifies  the  following  fundamental 
hypotheses: 

• Hypothesis  1 : The  spatiotemporal  dynamics  observed  during  gradient  sensing  are 
analogous  to  the  dynamics  of  the  activator-inhibitor  class  of  reaction-diffusion 
models. 

• Hypothesis  2:  The  membrane-resident  phosphoinositides  play  the  role  of  activator. 
We  can  address  both  these  hypotheses  by  experimentally  studying  the  features  of 

the  continuation  diagram  (Figure  2-5).  According  to  the  model,  when  the  mean 
chemoattractant  concentration  is  small  and  the  activator  synthesis  rate  is  slow,  the 
inhibitor  rapidly  counteracts  the  synthesis  of  the  activator  and  abrogates  its  polarization. 
When  the  chemoattractant  concentration  is  high  and  the  activator  synthesis  rate  is  fast; 
the  inhibitor  cannot  contain  the  synthesis  of  the  activator  and  causes  the  activator  to 
propagate  all  over  the  membrane.  Polarization  of  the  activator  can  occur  only  at 
intermediate  chemoattractant  concentrations,  when  the  synthesis  rate  of  the  activator  is 
optimum.  To  test  this  feature,  two  classes  of  experiments  have  been  performed  using 
neutrophil-like  HL60  cells  [72], 

• Biophysical  Experiments:  Cells  have  been  exposed  to  a range  of  uniform 
chemoattractant  concentrations  and  allowed  to  spontaneously  polarize.  The  extent 
of  polarity  observed  at  each  chemoattractant  concentration  has  been  quantified. 
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Biochemical  Experiments:  These  experiments  determine  whether  changes  in  PI5K 
activity  can  regulate  phosphoinositide  synthesis  and  influence  phosphoinositide 
polarization  in  a manner  consistent  with  the  model  and  with  our  choice  of  activator 
(membrane  phosphoinositides). 


Figure  4-1.  Kinetic  scheme  of  the  hypothesized  positive  feedback  loop. 

The  Positive  Feedback  Loop 

Positive  feedback  is  an  essential  component  of  the  gradient  sensing  model. 
Although  there  is  evidence  supporting  positive  feedback  [48],  the  components  of  this 
feedback  loop  are  yet  to  identified.  Based  on  the  current  literature,  one  can  construct  the 
model  shown  in  Figure  4-1. 

In  neutrophils,  the  receptors  for  the  chemoattractant  fMLP  are  coupled  to 
membrane-resident  G-proteins  (see  Figure  4-1).  When  a chemoattractant  molecule  binds 
to  a receptor,  the  G-protein  coupled  to  the  receptor  dissociates  into  its  Gar,  and  G Py 
subunits.  The  G Py  binds  and  activates  PI3K^  subunit,  resulting  in  the  synthesis  of 
PIP 3 from  PIP2  [42,  73],  The  PIP3  thus  produced  and  the  G Py  subunit,  then 


synergistically  activate  the  membrane-resident  protein,  P-Rexl,  which  belongs  to  the  Dbl 
family  of  Rac-GEFs  (guanine-nucleotide  exchange  factors  for  Rac)  that  activate  Rac  [45, 
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74],  In  resting  cells,  Rac  is  cytosolic  and  GDP-bound.  This  state  is  maintained  through  its 
interaction  with  a GDP  dissociation  inhibitor  (GDI)  (Figure  4-2).  Upon  receptor 
activation,  GDI  dissociates  from  Rac  by  an  unknown  mechanism,  and  Rac-GDP  targets 
to  membrane  via  its  C-terminal  prenylation,  where  Rac-GEFs,  such  as  P-Rexl,  activate 
Rac  by  dissociating  GDP  and  binding  GTP  [45],  To  be  sure,  other  Rac-GEFs  of  the  Dbl 
family  such  as  Tiaml  can  also  activate  Rac,  but  P-Rexl  is  most  important  insofar  as  its 
activity  represents  65%  of  the  total  Rac-GEF  activity  in  neutrophils  [45], 


Figure  4-2.  Components  of  the  positive  feedback  loop  in  a resting  cell  and  in  a polarized 
cell. 

In  vitro  experiments  have  shown  that  Rac  interacts  with  the  p85  subunit  of  PI3K  8 
activity  [46],  This  has  led  to  the  suggestion  that  PIP3  and  Rac  function  in  a positive 

feedback  loop  leading  to  rapid  accumulation  of  PIP3  [10,  48,  75],  However,  in  vivo  data 
show  that  less  than  1%  of  the  total  PI3K  activity  can  be  attributed  to  Rac-PI3K£  [76], 

A more  plausible  model  for  the  closure  of  the  positive  feedback  loop  is  the  Rac- 
mediated  activation  of  the  enzyme  PI5K,  which  catalyzes  synthesis  of  PIP  2 from  PI4P. 
Both  GDP-bound  and  GTP-bound  Rac  interact  with  PI5K  via  their  C-terminus  [47,  58, 
77,  78],  However,  only  GTP-bound  Rac  can  activate  PI5K[47,  78],  Based  on  these  facts. 
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the  following  model  is  proposed  (see  Figure  4-1,  Figure  4-2).  Under  resting  conditions,  a 
Rho-GDI/Rac-GDP/PI5K  complex  exists  in  the  cytosol.  Upon  receptor  activation,  this 
Rho-GDI  dissociates  from  the  complex  by  an  unknown  mechanism  and  the  Rac- 
GDP/PI5K  complex  gets  recruited  to  the  leading  edge  where  Rac-GEFs  convert  Rac- 
GDP/PI5K  to  Rac-GTP/PI5K.  Rac-GTP  may  either  activate  PI5K  or  provide  PI5K  access 
to  its  plasma  membrane  substrate,  PI4P.  This  results  in  localized  synthesis  of  PIP  2 at  the 

plasma  membrane.  This  increase  in  PIP2  drives  the  formation  of  more  PIP3  and  closes  the 
feedback  loop. 

The  forgoing  model  for  the  positive  feedback  loop  is  supported  by  several 
observations. 

• Introduction  of  membrane  permeant  PIP3  into  resting  neutrophils  triggers  the 
spontaneous  polarization  of  endogenous  PIP3  [48],  However,  if  PI3K  even 
membrane  permeant  PIP  3 cannot  induce  PIP3  polarization.  This  implies  that  a 

functional  positive  loop  is  required  for  polarization,  and  PI3K  is  an  essential 
component  of  this  loop. 

• Inhibition  of  PI3K  results  in  loss  of  directionality,  loss  of  PIP  3 localization  and 
improper  chemotaxis  of  neutrophils  [30,  75,  79], 

• Active  Rac-GTP  localizes  at  the  leading  edge  of  polarized  neutrophils[80]. 
Overexpression  of  Rac-GTP  results  in  PIP3  accumulation  all  over  the  plasma 
membrane  [49],  In  contrast,  Rac-GDP  inhibits  PIP3  accumulation  all  over  the 
plasma  membrane.  [49]. 

• In  human  fibroblasts,  PI5KIa  is  cytosolic  in  serum-starved  cells,  but  localizes  to 
membrane  ruffles  upon  growth  factor  stimulation  [58],  However,  inhibition  of 
either  PI5KIa  or  PI3K  completely  abolishes  ruffle  formation.  Overexpression  of 
Rac-GTP  and  PI5KI«  results  in  an  extensive  ruffling  phenotype,  whereas 
disruption  of  the  Rac-PI5KIa  interaction  results  in  inhibition  of  this  phenotype. 
Similar  effects  of  PI5K  on  actin  polymerization  have  been  observed  in  blood 
platelets  [47].  These  observations  are  consistent  with  our  model  where  PI3K,  PI5K 
and  Rac  are  essential  for  positive  feedback. 

• In  vitro  experiments  show  that  coexpression  of  PI5K  and  Rac-GTP  markedly 
stimulates  PI5K  activity,  whereas  coexpression  of  PI5K  and  Rac-GDP  inhibits 
PI5K  activity  [78,  81].  In  vivo  studies  show  that  PI5K  and  Rac-GTP  colocalize  at 
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the  plasma  membrane  and  induce  extensive  actin  polymerization.  In  contrast,  PI5K 
and  GDP-bound  Rac  fail  to  localize  at  the  membrane  and  completely  inhibit  actin 
polymerization  [81].  Thus,  Rac  regulates  both  the  activity  and  the  intracellular 
distribution  of  PI5K. 

• Recent  experiments  have  isolated  a novel  PI5K  homolog  called  PIPKH  [81].  This 
homolog  does  not  possess  PI5K  activity,  but  acts  as  a scaffolding  protein  that  binds 
to  PI5Ks.  Upon  overexpression,  PIPKH  results  in  small  increases  in  total  PIP2 
levels  and  dramatic  increases  in  PIP3  levels.  Consistent  with  the  idea  proposed  in 
our  model,  the  authors  suggest  that  PIPKH  recruits  PI5K  enzymes  to  specific 
intracellular  sites  where  small  increases  in  PIP2  levels  induces  massive 
accumulation  of  PIP3. 

• In  HL60  cells,  the  C-terminus  of  Racl  interacts  with  PI5K  [77],  Introduction  of  the 
purified,  membrane-permeant  Racl  C-terminus  (without  the  effector  region) 
inhibits  the  migration  of  HL60  cells.  Consistent  with  our  model,  this  inhibition 
occurs  because  the  endogenous  PI5K  is  saturated  with  inactive  Racl  C-terminus, 
and  therefore  cannot  synthesize  PIP2 . In  contrast,  when  a C-terminal  Racl  mutant 

that  does  not  bind  to  PI5K  is  introduced  into  HL60  cells,  motility  is  normal.  This 
observation  highlights  the  importance  of  the  Racl-PI5K  interaction  in  the  positive 
feedback  loop. 

Taken  together,  these  facts  strongly  suggest  the  positive  feedback  loop  has  the 
structure  shown  in  Figure  4-1. 

The  biochemical  experiments  described  in  this  chapter  check  if  variation  in  PI5K 
activity  regulates  phosphoinositide  polarity  in  a manner  consistent  with  the  model,  the 
hypothesized  positive  feedback  loop  and  the  choice  of  our  activator,  membrane 
phosphoinositides. 

HL60  cells,  a human  promyelocytic  leukemia  cell  line,  has  been  chosen  as  the 
model  system  to  perform  these  experiments.  When  treated  with  dimethyl  sulfoxide 
(DMSO)  for  6-7  days,  HL60  cells  differentiate  into  neutrophil-like  cells.  They  look  and 
behave  like  neutrophils,  orient  their  polarity  in  response  to  the  chemoattractant,  fMLP, 
and  migrate  in  the  right  direction  [18,  72],  The  differentiated  HL60  cells,  referred  to 
hereafter  as  cfHL60  cells,  are  easier  to  culture  than  the  terminally  differentiated 
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neutrophils.  Moreover,  HL60  cells  have  been  the  subject  of  numerous  gradient  sensing 
studies  providing  established  benchmarks  to  test  the  validity  of  protocols  used  in  this 
work. 

Protocols 

The  following  protocols  have  optimized  from  earlier  studies  [18,  22,  30,  48,  82], 
Culturing  and  Differentiation  of  HL60  cells 

HL-60  cells  (-O.lxlO6  #/ml)  are  grown  in  RPMI  1640  medium  (Sigma)  containing 
HEPES  buffer  and  supplemented  with  sodium  bicarbonate,  10%  FBS  (Atlanta 
Biologicals)  and  Penicillin-Streptomycin-Glutamine  (Invitrogen).  The  cells  are  allowed 
to  grow  in  8 ml  of  culture  medium  for  3-4  days  until  they  reached  a density  of  ~1.5xl  06 
#/ml.  Differentiation  of  the  HL-60  cells  to  neutrophils  is  then  induced  by  addition  of 
1.3%  DMSO  to  -3.0x1 06  cells  in  20  ml  culture  medium.  The  cells  are  propagated  for  7 
days  without  changing  the  medium.  Under  these  conditions,  the  suspended  cells 
transform  into  neutrophil-like  (/HL60  cells. 

Transient  Transfection  by  Electroporation 

An  electroporation  chamber  (Model  Series  1600,  BRL-GIBCO)  has  been  used  for 
transient  transfection.  Six  days  after  addition  of  DMSO,  -107  dHL60  cells  are  washed 
once  with  serum-free  RPMI-HEPES  and  resuspended  in  200p.l  of  the  same  medium.  For 
single  transfections,  30pg  of  plasmid  DNA  in  30pl  of  water  is  added  to  the  cells.  For 
double  transfections,  50pg  of  the  second  plasmid  is  included.  The  cell-DNA  mixture  is 
incubated  for  lOmins  at  room  temperature,  followed  by  5min  incubation  on  ice.  The 
mixture  is  transferred  to  the  electroporation  cuvettes,  and  then  subjected  to  an 
electroporation  pulse  at  330V,  lOOOpF  and  50Q.  The  transfected  cells  are  allowed  to 
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recover  at  room  temperature  for  10  minutes  and  subsequently  transferred  to  10  ml 
complete  culture  medium.  Transfection  efficiencies  varied  from  about  5%- 1 0% 
fluorescent  cells,  16-24  hrs  after  electroporation. 

Motility  Assays 
Preparation  of  cells 

Differentiated  HL-60  cells  are  washed  once  with  10  ml  of  serum  free  culture 
medium  and  resuspended  in  mHBSS  solution  (150mM  NaCl,  4mM  KC1,  1.2mM  MgCl2, 
10  mg/ml  glucose,  20mM  HEPES,  0.5%  low  endotoxin  human  serum  albumin,  pH  7.2, 
filter  sterilized).  The  cells  (~0.3  x 106  total)  are  plated  on  the  center  of  a sterile  no.  1 .5 
coverslip  rimmed  with  a square  agarose  spacer  (Labtek).  The  coverslip  was  previously 
coated  with  0.05mg/ml  human  fibronectin  (BD  Biosciences)  and  pre-blocked  with  0.5% 
human  serum  albumin.  After  30-40  minutes  of  incubation  in  a 37°C/5%  C02  incubator, 
the  non-adherent  cells  are  removed  by  two  washes  with  mHBSS/HSA.  Cells  are  then 
incubated  at  RT  for  5mins  in  200pl  of  mHBSS/HSA. 

Live  cell  imaging 

For  point  source  stimulation,  the  cells  are  exposed  to  a 0.5p.m  micropipette 
(Eppendorf  Femtotips)  backfilled  with  3pl  of  lOpM  fMLP  (Sigma)  in  mHBSS/HSA.  The 
micropipette  is  placed  in  the  microscope's  field  of  view  using  a joystick-controlled 
micromanipulator  (Zeiss).  If  necessary,  a small  backpressure  is  applied  to  the 
micropipette  using  a microinjector  (Eppendorf)  to  release  air  bubbles.  The  pipette  is  then 
placed  about  2-3  cell  lengths  away  from  a resting  cell  and  allowed  to  passively  diffuse 
chemoattractant.  Subsequent  events  were  captured  at  10s  intervals.  For  uniform 
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chemoattractant  stimulation,  20pl  of  lOnM  fMLP  in  mHBSS/HSA  is  gently  added  to 
coverslip  instead. 

Images  are  acquired  using  a cooled  charged  coupled  device  on  a multiwavelength 
multi-wavelength  wide-field  three  dimensional  Delta  Vision  Restoration  microscope 
(Olympus  1X70)  in  which  shutters,  filter  wheels,  focus  movement,  and  data  collection  are 
all  computer  driven.  Cells  are  imaged  using  a 60  x 1.4  N.A.  lens  (Olympus,  Lake 
Success,  NY)  and  n=1.514  immersion  oil.  For  fluorescence  imaging  of  living  cells,  GFP, 
CFP  or  YFP  signals  are  obtained  on  single  optical  sections  (0.25p.m)  near  the  bottom  of 
the  cells.  When  necessary,  differential  interference  contrast  (DIC)  images  of  the  cells  are 
simultaneously  captured  using  a Nomarski  system  optimally  aligned  for  our  microscope 
system.  Further  image  manipulation  is  done  using  Softworx,  IVE  or  Adobe  Photoshop 
Pro  software. 

Fixing  and  Immunostaining 

For  uniform  stimulation,  20pl  of  chemoattractant  (lOnM  fMLP  in  mHBSS/HSA)  is 
gently  added  to  the  coverslip.  Cells  are  allowed  to  respond  for  a stipulated  time  and  fixed 
for  20mins  with  3.7%  formaldehyde  in  freshly  prepared  cytoskeleton  buffer  containing 
lOmM  HEPES,  pH  7.2,  138mM  KC1,  3mM  MgCl2,  2mM  EGTA  and  0.32M  sucrose. 
Cells  are  then  washed  once  with  cytoskeleton  buffer  and  permeabilized  for  2 minutes 
using  0.2%  Triton-X  100  in  cytoskeleton  buffer  without  sucrose.  For  actin  staining,  cells 
are  incubated  for  one  hour  in  TxRed  conjugated  phalloidin  (Molecular  Probes)  in  PBS 
containing  1%  goat  serum,  washed  three  times  with  PBS  and  mounted  on  glass  slides. 
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Figure  4-3.  Spatiotemporal  dynamics  of  PIP3  localization  in  i/HL60  cells,  (a)  Image  of  a 
GFP-PH-Akt  transfected  c/HL60  cell  moving  towards  a micropipette  filled 
with  lOpM  sulphorhodamine  (red)  and  lOpM  fMLP.  The  GFP-PH-Akt 
localizes  to  the  leading  edge  of  the  migrating  cell,  (b)  A polarized,  GFP-PH- 
Akt  transfected  gML60  cell  migrating  in  response  to  a uniform  bolus  of 
lOOnM  fMLP.  (c)  A micropipette  filled  with  lOpM  fMLP  (white  circle)  is 
placed  towards  the  tail  of  a pre-pol ari zed  cell  (frame  1).  The  cell  responds  by 
turning  its  leading  edge  in  the  direction  of  the  micropipette  (frames  2,3,4)  ‘H’ 
denotes  the  head  of  the  cell.  ‘T’  denotes  the  tail. 

For  all  other  staining,  permeabilized  cells  are  incubated  in  a blocking  solution 
(blotto)  containing  50mM  Tris-HCl,  pH  7.5,  ImM  CaCb  and  3%  dry  milk,  followed  by 
one  hour  incubation  with  an  appropriate  dilution  of  the  primary  antibody  in  blotto.  After 
three  successive  washes  in  a solution  containing  25mM  Tris-HCl,  pH  7.4,  140mM  NaCl, 
2.5mM  KC1,  and  ImM  CaCh,  the  cells  are  incubated  for  one  hour  with  an  appropriate 
dilution  of  the  fluorescently  labeled  secondary  antibody.  After  three  successive  washes 
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with  the  wash  buffer,  the  stained  cells  are  mounted  on  glass  slides.  Fluorescence  images 


of  the  stained  cells  are  acquired  using  the  Delta  Vision  Restoration  microscope. 


Results 


Biophysical  Experiments 

Spatiotemporal  dynamics  of  gradient  sensing 

Previous  biophysical  studies  on  the  dynamics  of  the  PIP3  localization  in  c/HL60 
cells  have  shown  the  existence  of  spontaneous  polarization  and  polarized  sensitivity.  The 
validity  of  the  protocols  used  in  the  following  experiments  can  be  verified  by  first 
replicating  these  earlier  studies.  To  this  end,  HL60  cells  stably  transfected  with  a marker 
of  PIP3,  GFP-PH-Akt  (cell  line  from  UCSF,  Cell  culture  facility)  were  differentiated  with 
1 .3%  DMSO.  Seven  days  after  DMSO  treatment,  the  differentiated  cells  were  plated  on 
coverslips,  exposed  to  different  types  of  chemoattractant  profiles  and  the  resulting 
response  was  observed. 

• Response  to  a chemoattractant  gradient.  Cells  plated  on  the  coverslip  were 
exposed  to  micropipette  filled  with  lOpM  of  fMLP  in  mHBSS/HSA.  The  resting 
cell  polarized  in  morphology  and  accumulated  GFP-PH-Akt  at  its  leading  edge. 
This  was  followed  by  the  persistent  migration  of  the  cell  towards  the  pipette  tip 
(Figure  4-3a). 

• Response  to  uniform  chemoattractant  concentrations.  Cells  plated  on  the 
coverslip  were  exposed  to  a 100nM  uniform  concentration  of  fMLP  in 
mHBSS/HSA.  The  GFP-PH-Akt  construct  accumulated  at  the  plasma  membrane 
within  1 min.  Roughly  2 mins  later,  the  GFP-PH-Akt  polarized  spontaneously  at  a 
random  location  on  the  cell  membrane  (Figure  4-3b). 

• Response  to  a switch  in  the  direction  of  the  chemoattractant  gradient.  A 

micropipette  filled  with  lOpM  fMLP  was  placed  at  the  tail  of  a pre-polarized  cell. 
The  cell  exhibited  polarized  sensitivity  - its  leading  edge  turned  towards  the  pipette 
tip  (Figure  4-3c). 
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Spontaneous  polarization  as  a function  of  the  chemoattractant  concentration 

To  test  the  validity  of  the  continuation  diagram  (Figure  2-5),  c/HL60  cells  were 
starved  for  1 hour  in  mHBSS/HSA,  plated  on  coverslips  and  exposed  various  uniform 
fMLP  concentrations  varying  from  0.1  nM  to  lOOpM.  After  2 minutes  of  treatment,  the 
cells  were  fixed  with  3.7%  formalin,  permeabilized  and  stained  with  rhodamine 
phalloidin,  a marker  for  polymerized  actin. 


Figure  4-4.  Spontaneous  polarization  as  a function  of  chemoattractant  concentration. 
Cells  fixed  two  minutes  after  uniform  fMLP  treatment  and  stained  with 
rhodamine  phalloidin  to  visualize  actin  polymers.  Scale  bar,  20pm.  (a)  Low 
chemoattractant  (0.1  nM  fMLP).  Very  few  cells  showed  actin  polymerization 
at  the  membrane,  (b)  Intermediate  chemoattractant  (100  nM  fMLP).  Actin 
polymerization  is  confined  at  the  leading  edge  of  the  cell,  (c)  High 
chemoattractant  (lOOpM  fMLP).  Actin  polymerization  spreads  all  over  the 
periphery  of  the  cell. 

Consistent  with  model  predictions,  at  low  concentrations  of  chemoattractant 
(0.1  nM  fMLP),  cells  did  not  show  actin  polymerization  at  the  membrane  (Figure. 4-4a). 
However,  at  intermediate  concentrations  of  chemoattractant  (lOOnM  fMLP),  the  cells 
polymerize  actin  exclusively  at  their  leading  edges  (Figure  4-4b).  At  higher 
chemoattractant  concentrations  (lOOpM  fMLP),  cells  lost  their  tails  and  showed  actin 


polymerization  all  over  the  periphery  of  the  cell  (Figure  4-4c). 
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The  polarity  observed  in  each  of  these  experiments,  has  been  quantified  by 
measuring  the  ratio  of  the  cell  perimeter  showing  intense  phalloidin  staining  to  the  total 
cell  perimeter.  Results  of  this  quantification  have  been  shown  in  Figure  4-5. 


Figure  4-5.  Quantification  of  polarity.  ImageJ  software  was  used  to  the  measure  the  ratio 
between  the  cell  perimeter  showing  intense  phalloidin  staining  and  the  total 
cell  perimeter.  Results  are  from  two  independent  experiments  and  a total  of  30 
cells  were  counted  in  each  experiment. 

Biochemical  Experiments 

Role  of  PI5K  in  the  Positive  Feedback  Loop 

It  has  been  proposed  that  PI5K  is  an  element  that  completes  the  positive  feedback 
loop  for  PIP3and  PIP2  synthesis.  If  this  hypothesis  is  true,  inhibition  of  PI5K  activity 
should  block  the  localization  of  PIP3.  Overexpression  of  PI5K  activity  should  cause 
propagation  of  PIP3  all  over  the  membrane. 

Type  I PI5K  has  three  human  isoforms  PI5KIa,  PI5Kp  and  PI5Ky.  Based  on  the 
study  performed  by  [58],  the  PI5KIa  isoform  has  been  chosen  for  the  experiments  on 
HL60  cells.  This  study  showed  that  in  PDGF-stimulated  human  fibroblasts,  only  the 
PI5KIa  isoform  gets  recruited  to  the  membrane  and  plays  a role  in  membrane  ruffling. 
Both  the  other  isoforms  do  not  play  a role  in  actin  polymerization.  The  PI5K|3  isoform 
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localized  in  the  cytosolic  region  [58]  and  the  PI5Ky  isoform  associated  with  focal 
adhesions  [83]. 

Construction  of  fluorescent  PI5KIa  mutants  and  Coronin 

In  order  to  overexpress  or  inhibit  PI5KIa,  fluorescent  cyan  and  yellow  fusion 
proteins  of  full-length,  human,  PI5KIa  and  its  inactive  mutant,  PI5KIa£w  have  been 
constructed.  The  PI5KIafoz  mutant  has  two  mutations  in  its  catalytic  domain,  an  Arg-427 
to  Gin  (R427Q)  and  an  Asp-309  to  Asn(D309N).  This  inactive  mutant,  PI5KIa*w,  has 
been  shown  to  have  <1%  wild  type  PI5K  activity  [84],  PI5KIa£w  inhibits  PIP2  synthesis 
by  competing  with  the  endogenous  PI5KIoc,  for  the  substrate,  PI4P. 

The  full  length  PI5KIa  and  its  inactive  mutant  were  subcloned  by  PCR  from  the 
pcDNA3-pFLAG  vectors  (gift  from  Dr.  Richard  Anderson,  University  of  Wisconsin)  into 
pCI-Neo-eYFP  and  pCI-Neo-eCFP  mammalian  expression  vectors  using  BamHI  and 
EcoRI  cloning  sites.  All  the  PI5KIa  inserts  were  confirmed  by  sequencing.  The 
expression  of  YFP-PI5KIa  in  Cos7  cells  was  also  confirmed  by  immunoblotting  cell 
lysates  with  an  anti-GFP  antibody  (Chemicon).  In  order  to  observe  de  novo  actin 
polymerization  in  living  cells,  yellow  fusion  proteins  of  human  Coronin  1 have  also  been 
constructed.  Coronin  is  a 55  Kda  protein  that  was  isolated  from  the  slime  mould, 
Dictyostelium  [85],  Coronin  proteins  are  recruited  from  the  cytosol  to  a nascent  leading 
edge,  where  they  participate  in  actin  reorganization.  Appropriate  PCR  has  been  used  to 
insert  the  cDNA  for  the  most  abundant  human  Coronin  homologue,  Coronin  1 into  the 
pCI-Neo-eYFP  expression  system.  The  expression  of  YFP-Coronin  1 in  Cos7  cells  has 
also  confirmed  by  immunoblotting  with  the  anti-GFP  antibody. 
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Assaying  the  efficacy  of  the  fluorescent  constructs  in  Cos7  cells 


Figure  4-6.  Distribution  of  endogenous  PI5KIa  in  PDGF  stimulated  Cos7  cells,  (a) 

Endogenous  PI5KIa  colocalizes  with  (b)YFP-Coroninl,  at  membrane  ruffles 
(white  arrows). 
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Figure  4-7.  Ruffle  inhibition  in  PI5KIa£w  transfected  Cos7  cells.  Cells  expressing  both 
constructs,  (a)  CFP-PI5KIa^;  and  (b)  YFP-Coroninl,  showed  no  membrane 
ruffling  or  Coroninl  localization  In  contrast  cells  transfected  with  only  (c) 
YFP-Coroninl  showed  normal  membrane  ruffling  and  Coroninl  localization 
(white  arrows). 

Distribution  of  endogenous  PI5KIa  in  Cos7  cells.  Cos7  fibroblasts  were 
transfected  with  the  YFP-Coroninl  (marker  for  polymerized  actin)  construct  using  a 
Superfect  transfection  kit  (Qiagen).  These  cells  were  then  serum  starved  for  24  hours, 
stimulated  with  20  ng/ml  PDGF,  fixed  and  stained  with  a polyclonal  PI5KIa  antibody 
(gift  from  Dr.  Richard  Anderson,  University  of  Wisconsin).  Images  of  these  cells  show 
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that  endogenous  PI5KIa  (Figure  4-6a)  colocalized  with  YFP-Coroninl  (Figure  4-6b)  at 
membrane  ruffles. 


Figure  4-8.  Expression  and  distribution  of  PI5KIa  in  <iHL60  cells,  (a)  i/HL60  lysates 

were  probed  for  ubiquitous  expression  of  PI5KIa  using  1 pg  polyclonal  anti- 
PI5KIa.  A 68  kda  protein  corresponding  to  PI5KIa  was  found,  (b  and  c) 
GFP-PH-Akt  transfected  cell  dFIL60  cells  were  stimulated  with  lOOnM  fMLP. 
Cells  were  fixed  and  stained  with  3pg  of  anti-PI5KIa  antibody.  Both  GFP- 
PH-Akt  (b)  and  PI5KIa  (c)  enriched  ruffles  at  the  leading  edge  of  the  cell. 

PI5KIaArt  inhibits  ruffles.  [58]  showed  that  the  PI5KIai?r,  inhibits  actin 
polymerization  in  human  fibroblasts  treated  with  PDGF.  This  experiment  has  been 
reproduced  to  check  the  efficacy  of  the  PI5KIa£«  construct.  Cos7  cells  were 
cotransfected  with  YFP-Coroninl  and  the  inactive  mutant,  CFP-PI5KIaA«,  at  a 1:3  ratio. 
In  serum-starved  cells,  CFP-PI5KIaA?i  localized  predominantly  at  the  Golgi  (which  is 
consistent  with  the  results  in  [81]).  Upon  stimulation  with  20ng/ml  PDGF,  CFP-PI5KIa 
was  recruited  all  over  the  periphery  of  the  cell  (Figure  4-7a).  However,  87%  of  the  cells 
transfected  with  both  CFP-PI5KIo£h  and  YFP-Coroninl,  did  not  show  any  membrane 
ruffling  or  Coroninl  localization  (Figure  4-7a,b).  In  contrast,  71%  of  the  cells  transfected 
with  only  YFP-Coroninl,  showed  normal  ruffling  and  Coroninl  localization  (Figure  4- 
7c).  These  results  indicate  that  the  PI5KIafej  construct  is  functioning  properly. 
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Expression  and  distribution  of  PI5KIa  in  rfHL60  cells 

Before  studying  the  role  of  PI5KIa  in  dHL60  cells,  it  was  imperative  check  if  the  a 
isoform  of  Type  I PI5K  was  expressed  in  dHL60  cells.  Hence,  dHL60  cell  lysates  were 
probed  for  PI5KIa  using  the  polyclonal  PISKIa  antibody.  The  western  blot  isolated  a 
band  corresponding  to  the  68  kda  PI5KIa  (Figure  4-8a).  This  confirmed  the  expression  of 
PI5KIa  in  dHL60  cells. 

The  distribution  of  endogenous  PI5KIa  in  dHL60  cells  was  also  studied.  GFP-PH- 
Akt  transfected  stables  were  treated  with  lOOnM  chemoattractant  fMLP,  fixed  and 
stained  with  the  polyclonal  PI5KIa  antibody.  Both  GFP-PH-Akt  (Figure  4-8b)and 
PI5KIa  (Figure  4-8c)  were  predominantly  found  in  membrane  ruffles  at  the  leading  edge 
of  the  cell. 


Figure  4-9.  Overexpression  of  PI5K.  Cell  transfected  with  GFP-PH-Akt  and  CFP-PI5KIa 
was  exposed  to  lOOnM  uniform  fMLP.  (a)  GFP-PH-Akt  and  (b)  CFP-PI5KIa 
colocalized  all  over  cell  membrane,  (c)  Cell  transfected  with  only  GFP-PH- 
Akt,  localizes  the  probe  at  its  leading  edge.  Scale  bar,  20pm.  Graphs  in  the 
background  indicate  the  fluorescence  instensity  along  the  line  drawn  across 
the  cell. 
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PI5KIa  regulates  PIP3  polarity  in  </HL60  cells 

The  aim  of  this  experiment  was  to  check  if  regulation  of  PI5K  activity  could 
influence  PIP3  polarization.  GFP-PH-Akt  and  CFP-PI5KIa  (either  the  full  length  or  the 
inactive  construct)  were  cotransfected  into  r/HL60  cells  by  electroporation4,  six  days  after 
DMSO  treatment.  After  about  24  hours,  5-10%  of  the  cells  expressed  at  least  one  of  the 
constructs.  However,  nearly  75%  of  the  transfected  cells  expressed  both  the  constructs. 
Cells  were  then  plated  on  coverslips  and  stimulated  with  lOOnM  uniform  fMLP. 

Overexpression  of  PI5K.  When  cells  transfected  with  both  CFP-PI5KIa  and  GFP- 
PH-Akt  were  exposed  to  uniform  fMLP,  both  the  constructs  failed  to  polarize.  Instead, 
both  the  constructs  colocalized  all  over  the  periphery  of  the  cell  (Figure  4-9a,b).  These 
cells  were  neither  able  to  morphologically  polarize,  nor  migrate.  In  contrast,  cells  that 
were  transfected  with  only  GFP-PH-Akt,  polarized  normally  and  showed  GFP-PH-Akt 
localization  at  the  leading  edge  (Figure  4-9c). 


Figure  4-10.  Control  experiment.  Cell  transfected  with  GFP  (without  PH-Akt)  and  CFP- 
PI5KIa  was  exposed  to  lOOnM  uniform  fMLP.  (a)  GFP  (no  insert)  did  not 


4Direct  microinjection  of  the  constructs  into  the  cell  would  have  provided  a more  reliable  way  of 
transfecting  the  cells.  In  this  case,  expression  of  the  constructs  would  take  a short  time  and  would  have 
limited  side  effects  on  the  normal  function  of  the  cell.  However,  the  small  size  of  the  dHL60  cells  (~10pm) 
make  microinjection  difficult.  Hence,  electroporation  has  been  used. 
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localize  at  the  cell  membrane.  However, (b)  CFP-PI5KIa  colocalized  all  over 
cell  membrane.  Scale  bar,  11pm. 

The  colocalization  observed  in  Figure  4-9a,b  was  not  an  artifact  resulting  from 
plasma  membrane  folding.  This  was  confirmed  by  cotransfecting  dHL60  cells  with  CFP- 


PI5K-WT  and  GFP  (without  PH-Akt).  Although  CFP-PI5Ka  (Figure  4-1  Ob)  accumulated 


all  over  the  periphery  of  the  cell,  GFP  (Figure  4-10a)  remained  in  the  cytosolic  region. 


Figure  4-11.  Inhibition  of  PI5K.  Cell  transfected  with  GFP-PH-Akt  and  inactive,  CFP- 
PI5KIakn  was  exposed  to  lOOnM  uniform  fMLP.  (a)  GFP-PH-Akt  did  not 
localize  at  any  part  of  the  membrane  (b)  CFP-PI5KIa  localized  all  over  cell 
membrane  and  at  the  golgi.  Scale  bar,  20p.m.  Graphs  in  the  background 
indicate  the  fluorescence  instensity  along  the  line  drawn  across  the  cell. 
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Figure  4-12.  PI5K  regulation  of  PIP3  polarization  -Quantification  of  study. 
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Inhibition  of  PI5K.  Cells  cotransfected  with  the  inactive,  CFP-PI5KIakn  and 
GFP-PH-Akt,  were  exposed  to  fMLP  and  images  in  the  CFP  and  GFP  channels  were 
obtained.  CFP-PI5KIakn  localized  all  over  the  periphery  of  the  cell  and  at  in  the  golgi 
(Figure  4-1  lb).  However,  GFP-PH-Akt  failed  to  localize  at  any  portion  of  the  cell 
membrane  (Figure  4-1  la).  This  shows  that  PI5KIa  is  critical  for  the  polarization  of  GFP- 
PH-Akt.  Interestingly,  these  cells  were  still  able  to  polarize  and  migrate.  Cells  transfected 
with  only  GFP-PH-Akt,  showed  normal  polarization  of  the  probe  upon  treatment  with 
fMLP. 

Quantification  of  the  PI5K  studies  indicate  that  nearly  90%  of  the  cells  transfected 
only  with  GFP-PH-Akt  showed  normal  PIP3  polarization  (Figure  4- 12c).  However,  when 
PI5KIa  was  overexpressed,  nearly  80%  of  the  cotransfected  cells  localized  GFP-PH-Akt 
all  over  the  cell  membrane  (Figure  4- 12a).  In  contrast,  when  PI5KIa  was  inhibited,  65% 
of  the  cotransfected  cells  did  not  show  GFP-PH-Akt  localization  at  the  membrane  (Figure 
4- 12b). 

These  results  support  the  hypothesis  that  PI5KIa  is  involved  in  a positive  feedback 
loop  for  generation  of  PIP3.  They  also  indicate  that  PI5KIa  regulates  PIP3  polarization  in 
a manner  consistent  with  the  model  characteristics  shown  in  Figure  2-5. 

It  is  important  to  note  that  these  experiments  do  not  categorically  prove  the  role  of 
PI5KIa  in  a positive  feedback  loop  for  PIP3  synthesis.  This  hypothesis  can  be  verified  by 
inhibiting  PI5KIa,  and  then  introducing  membrane  permeant  PIP3  into  cells  to  check  if 
PIP3  can  still  trigger  its  own  synthesis. 


CHAPTER  5 
CONCLUSIONS 


A simple  reaction-diffusion  model  based  on  the  phosphoinositide  cycle  can 
effectively  capture  the  spatiotemporal  dynamics  involved  in  gradient  sensing.  The  model 
consists  of  two  subsystems:  (1)  Adaptation  subsystem  that  is  activated  on  a fast  time 
scale  and  inhibited  on  a slow  time  scale.  (2)  Polarization  subsystem  that  consists  of  a 
slow  diffusing,  autocatalytic  activator  and  a fast  diffusing  inhibitor.  Although  there  is 
some  evidence  supporting  the  choice  of  activator  (membrane  phosphoinositides)  and 
inhibitor  (inositol  phosphates),  more  experiments  need  to  be  performed  to  categorically 
prove  the  identity  of  these  species.  However,  the  success  of  this  model  suggests  that  the 
molecular  players,  whatever  be  their  identity  should  possess  the  properties  of  the 
variables  in  the  model.  The  features  exhibited  by  the  model  are  summarized  below. 

Fast  Responses 

1 . Phosphoinositide  localization.  Exposure  of  a cell  to  a shallow  chemoattractant 
gradient  results  in  the  translocation  of  phosphoinositides  to  the  leading  edge  of  the 
cell.  Receptor-ligand  binding  activates  the  model  and  nullifies  the  threshold 
required  for  the  amplification  of  membrane  phosphoinositdes.  This  results  in  the 
synthesis  of  phosphoinositdes  at  the  leading  edge.  The  rapid  diffusion  of  the 
inositol  phosphates  from  the  leading  edge  to  the  trailing  edge  inhibits  the  build  up 
of  phosphoinositides  at  the  back  of  the  cell  and  restricts  the  phosphoinositide 
localization  to  the  front  of  the  cell. 

2.  Unique  localization.  When  the  cell  is  exposed  to  two  unequal  sources  of 
chemoattractant,  phosphoinositides  localize  towards  both  sources  initially. 
However,  as  time  progresses,  only  the  localization  at  the  stronger  site  persists. 
Thus,  the  model  demonstrates  that  the  cell  can  migrate  only  along  a unique 
direction  regardless  of  the  chemoattractant  stimulus. 

3.  Properties  of  the  localization.  The  geometry  of  the  phosphoinositide  distribution 
varies  with  the  shape  of  the  chemoattractant  profile.  The  phosphoinositide 
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distribution  also  depends  on  the  mean  value  of  the  chemoattractant  concentration. 
Small  or  large  mean  values  yield  a nearly  uniform  phosphoinositide  distribution. 
Intermediate  mean  values  result  in  a pronounced  phosphoinositide  localization. 
The  width  of  the  phosphoinositide  peak  increases  with  the  mean  value. 

4.  Polarized  sensitivity.  When  a polarized  cell  is  exposed  to  a modest 

chemoattractant  gradient  along  a new  direction,  a new  localization  does  not 
develop  at  the  point  of  highest  chemoattractant  concentration.  Instead,  the  existing 
localization  turns  towards  the  new  stimulus.  If  the  new  gradient  is  localized  and 
steep  enough  to  overcome  the  existing  threshold,  a new  localization  is  generated 
and  the  existing  localization  is  disrupted.  If  the  polarized  cell  is  exposed  to  a new 
gradient  that  is  superimposed  on  a pre-existing  gradient,  the  response  observed 
depends  on  the  original  gradient,  the  intensity  of  the  pre-existing  polarity  and  the 
shape  of  the  new  gradient. 

Slow  Responses 

1 . Adaptation.  When  a cell  that  is  being  exposed  to  a uniform  chemoattractant 
concentration  is  subjected  to  a higher  uniform  chemoattractant  concentration, 
there  is  a rapid  increase  in  phosphoinositides  along  the  entire  cell  membrane.  This 
increase  is  followed  by  a slow  decline  in  the  phosphoinositides  back  to  its  initial 
distribution.  We  have  captured  this  response  by  assuming  that  a hypothetical 
component  gets  activated  on  a fast  time  scale  and  inhibited  on  a slow  time  scale. 
The  model  assumes  that  this  hypothetical  component  drives  phosphoinositide 
synthesis. 

2.  Spontaneous  polarization.  Experiments  show  that  cells  exposed  to  a uniform 
concentration  of  chemoattractant  polarize  before  adaptation  is  complete.  To 
account  for  this  property,  we  assume  that  the  chemoattractant  concentration 
undergoes  stochastic  fluctuations.  If  the  mean  concentration  of  the 
chemoattractant  is  large  enough  to  overcome  the  threshold,  the  uniform 
phosphoinositide  distribution  develops  into  a phosphoinositide  peak  at  a random 
location  on  the  cell  membrane.  As  time  progresses,  the  cell  adapts,  the  rate  of 
phosphoinositide  synthesis  decreases,  and  the  peak  disrupts. 

All  the  above  responses  clearly  preclude  the  hypothesis  that  the  cell  merely 
amplifies  the  external  signal.  Instead,  this  model  indicates  that  the  cell  must  be  viewed  as 
a system  that  non-linearly  processes  chemoattractant  inputs.  It  has  been  shown  in 
particular  that  these  complex  dynamics  can  be  explained  very  simply  in  terms  of  the 
instabilities  and  wavefront  dynamics  that  are  characteristic  of  the  activator-inhibitor 


model. 
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The  evolution  of  the  phosphoinositide  localization  consists  of  two  distinct  phases. 
During  the  first  phase,  the  region  of  the  membrane  that  crosses  the  threshold  gets 
amplified  to  a higher  steady  state.  All  other  regions  the  cell  get  inhibited  to  a lower 
steady  state.  The  formation  of  these  two  “plateaus”  result  in  a large  transition  layer  where 
there  is  accumulation  of  phosphoinositides.  This  accumulation  of  phosphoinositides 
drives  the  outward  motion  of  the  transition  layer  until  there  is  a zero  net-flux  of 
phosphoinositides  within  the  layer. 

A zeroth  order  approximation  to  the  steady  state  solution  has  been  be  derived  by 
utilizing  the  difference  in  diffusion  between  the  activator  and  inhibitor.  Perturbation 
techniques  were  used  to  arrive  at  separate  inner  and  outer  solutions  to  the  problem.  The 
two  solutions  were  then  patched  up  to  obtain  the  steady  state  response.  Geometrical 
properties  of  the  response,  such  as  height  and  width  have  also  been  derived. 

A Turing  stability  analysis  of  the  homogenous  resting  state  has  been  performed. 
This  analysis  gives  insight  into  the  conditions  required  for  the  homogenous  steady  to  lose 
stability  and  spontaneously  polarize  into  a phosphoinositide  localization.  It  has  been 
found  that  this  phenomenon  occurs  only  when  the  synthesis  rate  of  the  phosphoinositides 
is  at  an  optimum. 

Finally,  experiments  have  been  performed  to  justify  the  choice  of  this  class  of 
activator-inhibitor  models  to  explain  gradient  sensing.  This  has  been  done  by  testing  the 
following  model  result:  phosphoinositide  polarization  is  possible  only  when  the 
phosphoinositide  synthesis  rate  is  optimum.  This  result  has  been  tested  by  performing  (a) 
biophysical  experiments  where  c/HL60  cells  were  exposed  to  a range  of  uniform  fMLP 
concentrations,  and  fixed  after  2 minutes.  Actin  polarized  only  at  an  intermediate 
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concentration  of  chemoattractant  (lOOnM).  At  high  concentrations  (100pM),  actin 
polymers  were  observed  all  over  the  cell  membrane  and  at  low  concentrations  (0.1  nM) 
the  response  was  completely  lost,  (b)  biophysical  experiments  where  the  activity  of 
PI5KIa  in  c/HL60  cells  was  altered  to  study  the  direct  influence  of  phosphoinositide 
synthesis  on  phosphoinositide  polarization.  As  expected,  overexpression  of  PI5KIa 
resulted  in  the  propagation  of  phosphoinositides  all  over  the  plasma  membrane  and 
inhibition  resulted  in  abrogation  of  the  phosphoinositide  response.  This  study  also 
supports  the  involvement  of  PI5KIa  in  a positive  feedback  loop. 


APPENDIX  A 
NOMENCLATURE 


Table  A-l.  Iv 

odel  parameters. 

Parameter 

Definition 

ci 

Basal  synthesis  rate  of  inositol  phosphates 

cp 

Basal  synthesis  rate  of  phosphoinositides 

C 

Circumference  of  the  cell 

D, 

Lateral  diffusivity  of  inositol  phosphates 

Dp 

Lateral  diffusivity  of  phosphoinositides  in  plasma  membrane 

D, 

Lateral  diffusivity  of  phosphoinositides  in  endoplasmic  reticulum 

kf 

Rate  constant  for  receptor  mediated  phosphoinositide  formation 

kr 

Rate  constant  for  removal  of  phosphoinositides 

kp 

Rate  constant  for  basal  degradation  of  phosphoinositides 

ki 

Rate  constant  for  basal  degradation  of  inositol  phosphates 

k+ 

Rate  constant  for  receptor-ligand  association 

k~ 

Rate  constant  for  receptor-ligand  dissociation 

r 

Concentration  of  active  ligand  bound  receptor 

P 

Concentration  of  membrane  phsophoinositides 

i 

Concentration  of  inositol  phosphates 

Ps 

Concentration  of  phosphoinositides  in  endoplasmic  reticulum 

l 

Concentration  of  ligand 

rn,f 

Rate  of  formation  of  membrane  phosphoinositides 

rP,r 

Rate  of  removal  of  membrane  phosphoinositides 

'PS 

Degradation  rate  of  membrane  phosphoinositides 

ri,d 

Degradation  rate  of  inositol  phosphates 

R 

Radius  of  the  cell 

s 

Membrane  length  per  unit  cell  area 

a 

Time  scale  for  evolution  of  component  U 

b 

Time  scale  for  evolution  of  component  V 

r, 

Concentration  of  total  receptors 
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Table  A-2.  Greek  symbols 


Symbol 

Definition 

Vi 

Dimensionless  rate  of  basal  synthesis  of  inositol  phosphates 

Vp 

Dimensionless  rate  of  basal  synthesis  of  phosphoinositides 

Si 

Dimensionless  angular  diffusivity  of  inositol  phosphates 

5P 

Dimensionless  angular  diffusivity  of  phosphoinositides  in  plasma 
membrane 

K 

Dimensionless  angular  diffusivity  of  phosphoinositides  in  endoplasmic 
reticulum 

Kr 

Dimensionless  rate  constant  for  receptor  mediated  phosphoinositide 
formation 

kp 

Dimensionless  rate  constant  for  basal  degradation  of  phosphoinositides 

Ki 

Dimensionless  rate  constant  for  basal  degradation  of  inositol  phosphates 

K~ 

Dimensionless  rate  constant  for  receptor-ligand  dissociation 

P 

Dimensionless  concentration  of  active  receptor 

n 

Dimensionless  concentration  of  phsophoinositides  in  plasma  membrane 

i 

Dimensionless  concentration  of  inositol  phosphates 

n> 

Dimensionless  concentration  of  phosphoinositides  in  endoplasmic 
reticulum 

X 

Dimensionless  concentration  of  ligand 

T 

Dimensionless  time 

Tr 

Dimensionless  time  constant  for  receptor-ligand  binding 

Dimensionless  time  constant  for  evolution  of  U 

Ty 

Dimensionless  time  constant  for  evolution  of  V 

£ 

Angular  coordinate  scaled  between  0 and  1 

APPENDIX  B 

ORDERS  OF  MAGNITUDE  OF  PARAMETERS 

The  orders  of  magnitude  of  the  parameters  in  Table  2-1  were  based  on  the 

following  experimental  observations: 

1 . Under  receptor-activated  conditions,  the  turnover  time  of  membrane 
phosphoinositides  is  2 seconds  [51];  hence: 

k p 2 

kfp],krsp,  ~ lsec'1  =>  Kfp  - - 1.0 

kfsp, 


5.  Under  unstimulated  conditions,  the  turnover  time  of  membrane  phosphoinositides 
is  300  seconds  [51];  hence, 

cp/pt  ~ 0.0 lsec'1  ^ -0.01 

Ksp, 

6.  Under  unstimulated  conditions,  the  membrane  contains  10-30%  of  the 
phosphoinositides  in  the  endoplasmic  reticulum  (n  ~ 0. 1).  Since  the 
phosphoinositide  cycle  turns  slowly  under  these  conditions,  the  phosphoinositide 
concentration  is  determined  by  degradation  and  de  novo  synthesis.  It  follows  that, 

ICpK^Vp^Kp-QA 

7.  No  parameter  values  could  be  found  for  i/a  . It  was  assumed  that: 

V,=VP=  0.01 

This  is  tantamount  to  assuming  that  under  basal  conditions,  the  turnover  time  of 
the  inositol  phosphate  pool  is  on  the  order  of  100  seconds. 

8.  Under  basal  conditions,  the  inositol  phosphate  pool  constitutes  40%  of  the 
phosphoinositides  in  the  cell  (i  ~ 0. 1 ) [86];  hence: 

Kti  k , ~ 0. 1 

9.  The  lateral  diffusivity  of  cytosolic  1(1, 4, 5)  P,  is  2.8X10-6  cm2  sec"1  [86],  Since  the 
cell  diameter  is  on  the  order  of  10“3  cm  [87], 
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D 1C 2 

DJC 2 ~ 1 sec'1  =>  S , = ~ 1 .0 

krsp, 

10.  The  lateral  diffusivity  of  membrane  phospholipids  is  on  the  order  of  10~9  cm2 
sec'1  [15];  thus: 


Dp/C 2 ~ 10-3  sec"'  =>  Sp  = - 10-3 

1 1 . The  diffusivity  of  the  stored  phosphoinositides  has  been  assumed  to  be  the  same  as 
that  of  the  membrane  phosphoinositides: 


D„  !C 2 

Sa  =— p- 10'3 

* Ksp, 
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